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1 .  Introduction 

The  transition  from  laminar  to  turbulent  flow  in  boundary  layers  has  been  a 
subject  of  active  research  for  many  years.  The  intense  interest  in  this  subject 
stems  not  only  from  its  intrinsic  Importance  to  an  understanding  of  boundary- 
layer  behavior  but  also  from  important  technological  considerations  associated 
with  many  aerodynamic  and  hydrodynamic  flows,  and  the  need  to  predict  essential 
design  criteria  and  to  exert  some  measure  of  control  over  the  state  of  the 
boundary  layer.  Considerable  progress  has  been  made  in  revealing  the  fundamen¬ 
tals  and  evolutionary  nature  of  the  instability  and  transition  processes  that 
are  involved,  and  excellent  summaries  of  the  progress  made  to  date  can  be  found 
in  references  1-5.  However,  despite  the  progress  that  has  been  made,  the 
inherent  complexities  of  the  problem  have  prevented  the  prediction  of  transi¬ 
tional  behavior  with  a  degree  of  accuracy  required  for  practical  applications. 
This  situation  results  from  the  fact  that  there  are  many  factors  which  affect 
the  transition  process  and  much  more  needs  to  be  known  about  the  nature  of  the 
interaction  of  such  factors  as  freestreara  turbulence ,  surface  roughness,  surface 
curvature,  surface  temperature,  pressure  gradients,  Mach  number,  noise  and 
vibration  on  the  transition  process,  either  singly  or  in  combination.  All  of 
the  aforementioned  factors  have  been  investigated  to  some  degree  with  surface 
roughness  receiving  the  most  attention.  This  emphasis  on  surface  roughness 
stems  not  only  from  the  obvious  interest  in  maximizing  the  extent  of  laminar 
flow,  but  also  from  the  converse,  in  that  it  is  often  desirable,  for  example,  in 
wind  tunnel  tests  to  "artificially"  trip  the  boundary  layer  to  make  it  turbu¬ 
lent.  In  this  regard  it  is  desirable  to  "trip"  the  boundary  layer  as  effi¬ 
ciently  as  passible,  i.e.,  with  as  small  a  roughness  as  po33lble,  not  only  to 
avoid  contributing  to  the  form  drag  but  to  avoid  "distortion"  of  the  resulting 
turbulent  boundary  layer  (reference  6) ,  and  the  effect  on  test  data  resulting 
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from  such  "distortion. 11  Although  much  work  has  been  carried  out  on  roughness 
Induced  transition,  too  little  attention  has  been  given  to  characterizing  the 
state  of  the  resulting  boundary  layer. 

The  experimental  activity  dealing  with  the  effect  of  roughness  on  boundary 
layer  transition  at  both  incompressible  and  compressible  flow  speeds  spans  a 
period  of  approximately  forty  years,  and  may  be  considered  to  consist  of  two 
approaches:  one  in  which  the  activity  concerns  itself  with  the  location  of 
transition  and  its  dependence  on  parameters  characterizing  the  boundary  layer 
and  roughness,  and  the  other  concerned  with  details  of  the  flow  behavior  and  its 
relation  to  instability  processes  and  the  mechanism  of  transition.  Using  the 
former  approach,  considerable  data  was  obtained  for  both  two-dimensional  and 
three-dimensional  roughness  elements  in  which  correlations  of  Reynolds  numbers 
of  transition  with  roughness  size  and  position  and  boundary  layer  thickness  were 
obtained.  Excellent  accounts  of  this  approach  with  appropriate  bibliographies 
have  been  given  by  Smith  and  Clutter  (reference  7)  and  Tanl  (reference  8) . 
Although  the  parametric  approach  has  provided  the  basic  data  base  which  is  used 
for  present  day  criteria  involving  the  effect  of  roughness,  it  has  not  clarified 
the  basic  processes. 

In  general,  the  effect  of  roughness  is  to  induce  earlier  transition. 
However,  it  has  been  demonstrated  that  the  behavior  of  a  three-dimensional 
roughness  element  (references  8,  9)  is  markedly  different  than  that  for  two- 
dimensional  roughness  elements.  It  is  not  as  effective  as  a  two-dimensional 
roughness  in  inducing  transition,  and  its  behavior  is  much  more  critical. 
Transition  moves  forward  gradually  with  increasing  velocity  for  a  two- 
dimensional  roughness,  whereas  for  a  three-dimensional  roughness  transition 
moves  rapidly  toward  the  roughness  with  a  relatively  small  increase  in  velocity 
after  a  critical  velocity  has  been  reached.  It  is  generally  recognized  that  in 
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order  to  materially  advance  the  state-of-the-art  and  our  understanding  of  the 
basic  mechanisms  Involved  in  this  behavior,  the  approach  of  investigating  the 
steady  and  fluctuating  flow  fields  in  detail  is  required.  This  more  basic 
approach  has  been  used  with  some  success  in  connection  with  two-dimensional 
roughness  elements,  and  it  has  been  shown  (reference  10)  that  the  basic 
mechanism  by  which  a  two-dimensional  roughness  element  induces  earlier 
transition  to  turbulent  flow  is  by  the  destabilizing  Influence  of  the  flow  In 
the  immediate  downstream  vicinity  of  the  roughness.  The  flow  downstream  of 
isolated  three-dimensional  roughness  elements  has  received  considerable 
attention  (references  11-21).  However,  much  of  this  work  utilized  flow 
visualization  techniques  and  the  resulting  insight  into  the  flow  processes  is 
primarily  qualitative  In  nature.  Detailed  measurements  of  the  flow  downstream 
of  a  three-dimensional  roughness  are  limited  and  as  yet  no  clear  explanation  of 
its  behavior  in  inducing  transition,  similar  to  that  for  two-dimensional 
roughness,  exists.  Nevertheless,  it  is  felt  that  the  behavior  should  in  some 
manner  be  stability  governed.  It  was  with  the  objective  of  providing 
Information  on  this  aspect,  to  considerably  extend  the  data  base,  and 
characterize  the  state  of  the  boundary  layer,  if  only  for  a  single  roughness 
element,  that  the  present  Investigation  was  undertaken. 

2.  Experimental  Arrangement  and  Procedure 
2.1  The  Flow  Facility 

The  Investigation  was  conducted  in  the  National  Bureau  of  Standards 
Unsteady  Flow  Facility  shown  schematically  in  figure  1.  The  facility  is  of  the 
closed-return  type  with  a  square  test  section,  1.37  m  in  cross-section  and 
4.88  m  in  length.  It  Is  powered  by  a  100  h.p.  constant  speed  motor  and  single 
stage  axial  fan  equipped  with  shutters  for  varying  the  flowrate.  The  tunnel  can 


9 


AEDC-TR-87-7 


be  operated  either  in  a  steady  mode  as  a  conventional  wind  tunnel,  or  in  an 
unsteady  mode  for  generating  oscillatory  flows  with  varying  frequency  and  ampli¬ 
tude.  Both  modes  of  operation  were  used  In  the  present  investigation. 

The  unsteady  operation  of  the  wind  tunnel  is  accomplished  by  using  a  system 
of  rotating  shutters  and  a  secondary  or  bypass  test  section.  The  secondary  test 
section  is  solely  for  the  purpose  of  providing  a  nearly  constant  volume  flowrate 
through  the  fan.  This  is  achieved  by  balancing  the  mean  flow  between  the  two 
sections  with  the  trim  shutters  and  by  alternatively  directing  part  of  the  flow 
between  the  two  sections  by  rotating  the  shutters  in  each  section  90°  out-of- 
phase.  Guide  vanes  are  installed  at  each  corner,  and  the  solid  lines  at  the 
guide  vane  positions  shown  in  the  figures  serve  to  indicate  the  beginning  and 
ending  of  the  partitioning  of  the  flow  in  the  wind  tunnel  circuit.  A  screen 
from  position  C  and  another  from  position  D,  as  indicated,  and  extending  to  the 
leading  edge,  AB,  of  the  dividing  partition,  assist  in  diverting  the  flow  into 
the  two  ducts  as  required  by  the  rotating  shutters.  Each  bank  of  shutters 
consists  of  nine  counter-rotating  shutters,  any  one  of  which  may  be  positioned 
at  a  fixed  angle  between  0°  and  90°.  By  varying  the  number  of  rotating 
shutters,  the  angle  of  the  fixed  shutters,  and  the  flowrate  through  the  fan, 
both  the  mean  velocity  and  the  velocity  amplitude  in  the  working  test  section 
can  be  controlled.  The  rotational  speed  of  the  shutters  determines  the 
frequency  which  ranges  from.O.T  Hz  to  25  Hz. 

The  lower  duct,  resulting  from  the  partitioning  of  the  flow,  is  the  primary 
duct.  It  contains  the  working  test  section  and  conforms  to  conventional  wind 
tunnel  design.  The  contraction  ratio  is  5.3,  and  the  settling  chamber  contains 
a  honeycomb  and  screens  for  flow  conditioning,  i.e.,  for  obtaining  a  spatially 
uniform  mean  flow,  and  a  reduced  freestream  turbulence  in  the  teat  section.  The 
maximum  velocity  that  can  be  attained  in  the  test  section  when  operating  in  a 
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conventional  steady  mode,  with  the  upper  duct  closed,  la  27  m/a.  The  centerline 
turbulence  intensity,  i.e.,  the  ratio  of  the  rms  value  of  the  axial  component  of 
the  velocity  fluctuations,  u^,  to  Uq,  a  reference  mean  velocity,  and  its 
variation  with  UQ  at  0.5  m  upstream  from  the  boundary  layer  plate  is  shown  in 
figure  2.  Despite  the  complexity  of  the  wind  tunnel  configuration,  the  tur¬ 
bulence  intensity,  above  12  m/s,  is  about  0.06  percent  and  is  characteristic  of 
turbulence  levels  of  conventional  low  turbulence  wind  tunnels.  Below  12  m/s, 
and  over  the  speed  range  covered  with  the  roughness  element  the  turbulence 
Intensity  increases  with  decreasing  velocity  to  about  0,15  percent  at  5  m/s. 

2.2  The  Flat  Plate  and  Traversing  Arrangement 
The  investigation  was  carried  out  using  a  boundary  layer  on  an  aluminum 
plate,  3.66  m  long,  1.22  m  wide  and  9.5  mm  thick  which  had  a  leading  edge  in  the 
form  of  a  half-wedge  with  a  5®  angle  on  the  non-working  side  of  the  plate.  The 
leading  edge  itself  is  rounded  with  a  radius  of  approximately  0.H  mm.  The  plate 
is  mounted  vertically  and  centrally  with  the  leading  edge  at  a  distance  1.0  m 
downstream  of  the  entrance  to  the  test  section.  A  false  wall  made  of.  plexiglass 
mounted  on  the  tunnel  wall  opposite  the  working  side  of  the  plate  permitted  the 
adjustment  of  the  pressure  gradient  within  moderate  limits  and  was  used  to 
compensate  for  the  boundary  layer  growth  and  establish  a  zero  pressure  gradient 
along  the  plate.  A  positive  angle  of  attack  was  obtained  by  some  blocking  of 
the  air  passage  on  the  working  side  of  the  plate  with  screens  of  1 6-mesh  and 
wire  diameter  0.28  mm.  Since  the  blocking  screens  had  to  be  placed  farther 
downstream  than  any  usable  portion  of  the  plate,  the  passage  was  extended  0.86  m 
by  joining  the  trailing  edge  of  the  plate  with  a  section  of  plywood.  This 
extension  also  prevented  a  sharp  change  in  pressure  along  the  plate  that  would 


11 


AEOC-TR-87-7 


otherwise  occur  once  the  carriage  for  the  traversing  equipment  passed  the  down¬ 
stream  end  of  the  plate.  The  positive  angle  of  attack,  with  the  stagnation 
point  on  the  working  side  of  the  plate,  provided  a  smooth  flow  condition  at  the 
leading  edge,  with  a  resulting  favorable  pressure  gradient  over  the  first  2  feet 
of  surface.  The  distribution  of  mean  velocity,  U1 ,  In  the  free-stream  outside 
the  boundary  layer,  with  distance  from  the  leading  edge,  x,  obtained  by  adjust¬ 
ing  the  false  wall,  is  shown  in  figure  3-  It  is  seen  that  beyond  the  first  two 
feet  of  surface  a  reasonable  zero  pressure  gradient  condition  has  been  obtained. 

A  view  of  the  experimental  arrangement  from  upstream  of  the  plate  and 
looking  downstream  Is  shown  In  figure  4.  As  can  be  seen,  the  surface  is  highly 
polished  to  a  mirror-like  finish.  Also  shown  in  the  figure  is  an  overall  view 
of  the  arrangement  for  traversing  the  boundary  layer.  Close-up  views  of  the 
roughness  element,  hot-wire  probe  and  traversing  mechanism  used  to  traverse 
either  in  a  spanwlse  direction  or  in  a  direction  normal  to  the  plate  surface  are 
shown  in  figures  5a  and  5b.  The  traverse  is  mounted  to  a  streamlined  strut 
which  can  ride  along  the  plate  on  two  rails  attached  at  the  top  and  bottom  of 
the  plate  for  traversing  in  a  downstream  direction.  The  measuring  probe,  which 
extends  25  cm  upstream  from  the  point  of  support  to  avoid  interference  effects, 
is  driven  by  a  micrometer  screw  of  0.5  mm  pitch.  Additional  probes  can  be 
attached  at  fixed  locations,  as  for  example,  the  Pitot-static  tube  in  figure  5a, 
for  measuring  the  free-3tream  velocity  during  boundary  layer  surveys.  The 
initial  distance  of  the  probe  from  the  surface  is  obtained  by  making  use  of  the 
mirror-like  finish,  and  using  a  prism  to  reflect  the  probe  and  its  image  onto  a 
calibrated  reticule  of  a  microscope. 


12 


AEDC-TR-8  7-7 


2.3  The  Laminar  Boundary  Layer 

It  was  necessary  before  proceeding  with  the  study  of  the  behavior  of  a 
roughness  element  to  determine  the  extent  and  character  of  the  laminar  boundary 
layer  that  was  established  on  the  plate.  The  extent  of  the  laminar  boundary 
layer  "window”  is  shown  in  figure  6.  The  figure  shows  the  well-known  phenomenon 
of  transverse  contamination  of  turbulent  flow  originating  at  the  top  and  bottom 
corners  of  the  plate  at  the  leading  edge.  The  transverse  contamination  angle, 
as  measured  from  the  rail  support  for  the  strut  is  about  10°,  in  good  agreement 
with  existing  data  for  such  behavior.  Figure  7  shows  the  mean  velocity  dis¬ 
tributions  across  the  boundary  layer  at  the  centerline  of  the  plate  at  two  free- 
stream  speeds,  and  for  two  downstream  positions.  The  two  downstream  positions 
were  chosen  so  as  to  cover  a  considerable  extent  of  the  laminar  boundary  layer 
"window."  The  distributions  are  plotted  nondimensionally  for  comparison  with 
the  Blaslus  distribution  for  a  zero-pressure  gradient  boundary  layer,  where  U  is 
the  x-component  of  the  local  mean  velocity  in  the  boundary  layer,  y  is  the 
distance  from  the  surface,  and  v  is  the  kinematic  viscosity.  It  may  be  that  the 
tendency  of  the  velocity  profiles  to  be  slightly  steeper  in  shape  than  the 
Blaslus  reflects  the  effect  of  the  favorable  pressure  gradient  over  the  forward 
part  of  the  plate,  and/or  a  virtual  origin  other  than  at  the  leading  edge,  but 
the  effect,  if  any,  is  sufficiently  small,  and  within  the  experimental  uncer¬ 
tainty  the  boundary  layer  may  be  considered  to  be  of  the  Blasius  type.  Conse¬ 
quently,  In  characterizing  the  flow  or  flow  parameters  at  the  roughness  element, 
with  the  roughness  element  present,  the  Blasius  distribution  is  used.  Although 

this  assumption  neglects  the  distortion  in  the  Immediate  vicinity  upstream  of 

1  /2 

the  roughness,  and  at  the  roughness,  the  Blasius  length  scale,  n  *=  (v  x/U^  , 
13  a  realizable  length  scale,  and  provides  consistency  with  other 
investigations.  Measurements  of  the  disturbance  level  in  the  boundary  layer 
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without  roughness  within  the  laminar  boundary  layer  "window"  are  shown  In  figure 
8.  Distributions  across  the  boundary  layer  of  the  rms  value  of  the  longitudinal 
component  of  the  velocity  fluctuation,  u' ,  are  shown  relative  to  its  value  in 
the  freestream,  uj ,  for  the  same  x-positions  and  for  about  the  same  freestream 
velocities  as  for  figure  7-  The  nature  of  the  distributions  is  such  that  one 
can  infer  that,  to  a  considerable  degree,  the  disturbance  level  is  a  consequence 
of  variations  of  boundary  layer  thickness  in  space  and  time  associated  with 
free-stream  turbulence  effects  as  proposed  in  reference  22  and  not  as  a  result 
of  a  Tollmien-Schlichting  type  of  instability.  The  maximum  in  the  disturbance 
level,  at  about  0.5  the  boundary  layer  thickness,  which  in  Blasius  coordinates 
is  at  2.5,  and  the  relatively  large  low  frequency  content  of  the  disturbance, 
i.e.,  below  the  Tollmlen-Schli cting  range  of  instability,  are  consistent  with 
this  conclusion.  This  behavior  is  illustrated  in  figure  9  which  shows  the 
spectral  content  relative  to  the  range  of  Tollmien-Schlichting  frequencies  for 
the  highest  boundary  layer  Reynolds  number  of  figure  8,  i.e.,  at  -  12.1  m/s 
and  x  -  21 3.H  cm,  for  which  the  boundary  layer  Reynolds  number  based  on 

i 

displacement  thickness,  ,  was  2070.  The  spectral  distribution  was. obtained 
in  the  region  of  the  maximum  in  u’/u*  at  n  *=  2.45,  and  the  Roman  numerals  I  and 
II  noted  on  the  figure  indicate  the  frequencies  corresponding  to  Branches  I  and 
II  of  the  stability  diagram  for  a  Blasius  flow.  The  nominal  band-width  of  the 
spectral  analyzer  is  proportional  to  the  band-pass,  and  for  the  100  Hz  band-pass 
it  is  0.3  Hz.  In  addition,  the  nominal  low-frequency  cut-off  of  the  analyzer  is 
0.4  Hz.  Consequently,  no  particular  significance  should  be  attached  to  the 
fall-off  in  spectral  intensity  that  is  shown  in  the  range  from  0.0  to  1.0  Hz. 

The  peak  in  the  spectral  Intensity  at  20  Hz,  in  view  of  its  proximity  to  Branch 
I,  should  not  be  Interpreted  as  being  due  to  a  Tollmien-Schlichting  instability. 
It  is  a  consequence  of  a  vibration  induced  by  the  fan  motor,  and  is  accentuated 
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by  the  velocity  gradient  in  the  boundary  layer.  Attempts  to  eliminate  the  20  Hz 
vibration  by  soft-mounting,  damping,  etc.  were  not  successful.  Although  it  was 
minimized  to  the  extent  shown,  it  was  not  possible  to  eliminate  it  inasmuch  as 
the  20  Hz  vibration  from  the  fan  motor  was  transmitted  to  the  flat  plate  and 
wind  tunnel  walls  as  evidenced  by  accelerometer  measurements.  The  free-stream 
disturbance  level,  u^ /U^ ,  along  the  plate  is  about  a  factor  of  1.9  ±  10  percent 
larger  than  the  corresponding  disturbance  level,  uVUo,  upstream  of  the  plate. 
This  increase  in  disturbance  level  along  the  plate,  although  it  has  never  been 
thoroughly  evaluated,  has  been  observed  in  previous  investigations.  However,  It 
is  not  unexpected,  and  may  be  attributed  to  potential  fluctuations  arising  from 
the  turbulent  boundary  layers  on  the  walls  of  the  tunnel.  The  magnitude  of  the 
effect  may  therefore  differ  from  one  flow  configuration  to  another.  In  addition 
there  may  be  very  low  frequency  potential  fluctuations,  i.e.,  below  2  Hz,  in  the 
free-stream  arising  from  an  intermittent  separation  in  the  wind  tunnel  circuit. 
It  Is  believed  that  this  behavior  may  be  responsible  for  much  of  the  Increase  in 
disturbance  Intensity  at  the  lower  velocities  shown  in  figure  2.  Another 
contributing  component  to  the  free-stream  disturbance  level  is  the  free-stream 
turbulence,  i.e.,  turbulent  vorticity  fluctuations,  which  were  the  type  of 
disturbances  considered  in  reference  22. 

It  is  evident  that  sound  and  vibration  are  factors  that  may  affect  the 
state  of  the  boundary  layer. and  are  important  considerations  particularly  in 
investigations  pertaining  to  flow  Instability.  In  order  to  assess  this  aspect 
the  spectrum  of  the  u-fluctuatlon  was  measured  in  the  free  stream  at  x  ■  221  cm, 
y  =  3.49  cm  and  =8.1  m/s.  The  spectrum  in  terms  of  the  relative  rms 
amplitude  is  3hown  in  figure  10,  No  particular  significance  should  be  attached 
to  the  relative  amplitudes.  The  position  at  which  the  spectrum  was  measured  was 
arbitrarily  selected,  and  the  amplitude  of  the  various  frequencies  may  vary  with 
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position  in  the  teat  section.  The  purpose  is  to  illustrate  the  spectral  peaks 
present  in  the  free-strean  turbulence  that  may  serve  as  a  source  of  input 
disturbances  for  any  flow  Instability.  The  source  of  the  spectral  peaks  is 
attributed  to  the  wind  tunnel  fan.  The  peak  at  20  Hz,  as  alluded  to  previously, 
corresponds  to  the  fan  motor  rpm.  The  peaks  at  180  Hz  and  the  higher  harmonics 
of  1 80  Hz  are  a  consequence  of  the  passage  frequency  of  a  9-bladed  fan.  A  more 
direct  evaluation  is  given  by  figure  11  which  shows  the  spectrum  of  the  output 
of  an  accelerometer  attached  to  the  fan  housing.  It  is  seen  that  the  peaks 
correlate  well  with  the  peaks  observed  in  the  spectrum  of  the  Tree-stream 
turbulence.  Microphone  spectra  taken  in  the  room  near  the  test  section  also 
corroborated  that  the  spectral  peaks  were  a  consequence  of  sound  and  vibration. 

It  is  noted,  as  is  reflected  in  figure  8,  that  there  is  an  anomaly  in  the 
increase  in  u'/u|  from  x  =  91.4  cm  to  x  -  213-4  cm  for  the  lower  free-stream 
velocity  as  compared  to  that  at  the  higher  free-stream  velocity.  There  is  no 
ready  explanation  for  this  behavior.  One  may  speculate  that  In  addition  to  the 
complex  interactions  referred  to  above,  the  situation  is  further  complicated  by 
the  transverse  contamination,  and  that  the  effect  of  free-stream  turbulence  on 
the  boundary  layer  may  also  depend  on  the  thickness  of  the  layer  relative  to  the 
scale  of  the  turbulence.  Be  that  as  it  may,  the  state  of  the  layer  is  such  as 
to  not  prevent,  nor  unduly  complicate,  an  adequate  study  of  the  effect  of 
roughness . 

2.4  Procedure 

The  roughness  elements  were,  for  the  mo3t  part,  hemispherical  with  nominal 
diameters,  d,  of  3.18  and  6.36  mm.  They  were  individually  studied  and  were 
attached,  using  rubber  cement,  to  the  surface  on  the  centerline  (z  =■  0)  at 
varying  distances  from  the  leading  edge,  x^,  which  were  well  within  the  laminar 
boundary  layer  "window”.  The  actual  roughness  height,  k,  due  to  the  method  of 
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attachment  and  the  tolerance  in  the  nominal  values,  could  differ  from  the 
nominal  value  and  was  measured  in  each  case  with  a  dial  indicator.  A 
cylindrical  element  with  d  =  3.18  mm,  k/d  =  1.0,  and  positioned  at  x^  -  91.4  cm, 
with  its  axis  in  the  y-direction  was  also  used. 

In  documenting  the  behavior  of  a  single  roughness  element  a  number  of 
different  aspects  required  attention.  One  was  the  variation  of  the  Reynolds 
number  of  transition,  Ret,  with  roughness  Reynolds  number,  Re^,  defined  by 
U  x  /v  and  U  k/v  respectively  where  x.  is  the  distance  from  the  leading  edge  of 

1  t  K  Zr 

the  plate  to  the  beginning  of  transition  to  turbulent  flow  in  the  boundary 
layer,  and  Uk  is  the  mean  velocity  that  would  exist  in  the  boundary  layer  at  the 
height  of  the  roughness  without  the  roughness  present.  Another  aspect  is  the 
parametric  behavior  of  the  frequency,  f,  of  eddies  generated  by  the  roughness. 
Attention  was  also  given  to  the  nature  of  the  boundary  layer  distributions  of 
mean  velocity  and  intensity  of  disturbances  at  various  distances  downstream  from 
the  roughness,  x,  as  well  as  to  the  effect  of  an  oscillatory  flow  on  eddy  shed¬ 
ding  and  roughness  Induced  transition.  In  regard  to  the  latter  two  aspects,  the 
behavior  of  only  one  size  of  roughness  at  one  roughness  position  was  investi¬ 
gated,  l.e.,  d  =  3-18  mm,  k  =  1.7  mm  and  =  91.4  cm. 

The  distributions  of  U  and  u' ,  at  x  ranging  from  1.27  cm  to  30.5  cm,  were 
measured  on  the  centerline  and  in  the  spanwise  direction  for  various  unit 
Reynolds  numbers,  using  conventional  constant  temperature  hot-wire  anerao- 

metry.  During  the  boundary  layer  surveys,  a  given  unit  Reynolds  number  was 
constant  to  within  ±  two  percent.  Platinum  wires  1.27  um  In  diameter  and  0.25 
mm  in  length  were  used,  and  no  correction  for  wire  length,  nor  for  the  nonlinear 
response  of  the  hot-wire  was  deemed  necessary.  The  length  of  the  wire  was 
selected  as  being  adequate  for  the  size  of  the  roughness  under  study  from 
comparative  measurements  of  the  mean  velocity  in  the  immediate  downstream 
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vicinity  of  the  roughness  with  different  wire  lengths.  Equipment  cufltomarily 
used  in  association  with  hot-wire  anemoraetry  for  processing  data  such  as  a  real¬ 
time  spectral  analyzer  with  an  on-line  x-y  plotter,  film  recordings  from  an 
oscilloscope,  and  strip-chart  recordings  also  proved  useful. 

3.  Results  and  Discussion 
3.1  Critical  Reynolds  Number 

The  existence  of  a  critical  roughness  Reynolds  number  for  correlating  the 
Reynolds  number  of  transition  for  three-dimensional  roughness  elements  was  first 
prepared  by  Schiller  (reference  23),  It  was  predicated  on  the  analogy  to  the 
critical  Reynolds  number  for  an  obstacle  in  uniform  flow  at  which  the  downstream 
flow  becomes  altered  by  vortex  shedding.  Although  considerable  data  on  the 
critical  roughness  Reynolds  number  for  three-dimensional  elements  exists  in  the 
literature  for  a  variety  of  geometrical  shapes  such  as  spheres,  cones,  and 
cylinders  (references  7,9,11,13,24),  no  comparable  data  exists  for  hemispherical 
roughness. 

The  critical  behavior  of  three-dimensional  hemispherical  roughness  elements 
in  inducing  earlier  transition  to  turbulent  flow,  as  observed  in  the  present 
study,  is  illustrated  in  figure  12.  The  measurements  shown  were  made  for  two 
different  sizes  of  roughness  mounted  at  different  positions  from  the  leading 
edge.  The  position,  xfc,  was  determined  by  placing  a  hot-wire  probe  close  to  the 
surface  along  the  centerline  at  different  locations,  x,  and  slowly  increasing 
until  turbulent  bursts  first  occurred.  The  results  clearly  demonstrate  the 
critical  behavior  of  the  transition  process  induced  by  the  roughness  in  that 
transition  moves  rapidly  upstream  toward  the  roughness  element  with  a  relatively 
small  increase  in  velocity  after  a  critical  velocity  has  been  reached.  The 
limiting  case  of  xfc  =  x^  is  shown  by  the  solid  line,  and  the  critical  Rek  is 
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determined,  as  in  reference  13,  by  the  extrapolated  intersection  with  the  limit¬ 
ing  curve  indicated  by  the  dashed  line.  The  Increase  in  Re^  at  the  lower  values 
of  Ret,  indicating  that  the  transition  is  not  moving  forward  as  rapidly,  is 
generally  observed  to  be  the  case  when  is  less  than  30-5  cm  from  the  rough¬ 
ness  element.  This  behavior  is  apparently  a  consequence  of  the  finite  distance 
required  for  the  transition  to  turbulence  associated  with  the  stability  of  the 
flow.  It  illustrates  that  considerable  care  should  be  given  to  the  state  of  the 
"tripped"  boundary  layer,  particularly  when  using  "trips"  in  wind  tunnel  test¬ 
ing,  in  order  to  ascertain  the  degree  to  which  the  boundary  layer  Is  free  of 
"distortion”.  This  question  is  addressed  further  for  a  given  size  of  roughness 
in  Section  3*4. 

It  should  be  noted  that  Re^  is  constant  with  a  value  of  about  325  for  the 
various  conditions  shown  in  figure  12.  This  is  at  variance  with  the  results  of 
reference  13  for  cylindrical  roughness  elements  for  which  the  critical  value  of 
Re^  was  found  to  vary  from  600  to  1000  and  Increasing  with  decreasing  Re  . 
However,  in  reference  13,  x^  was  defined  as  the  position  where  the  turbulent 
"bursts”  had  an  intermittency  factor  of  50  percent,  and  the  aspect  ratio,  k/d, 
of  the  cylindrical  elements  was  1.0  compared  to  a  value  of  0.5  for  hemispherical 
elements.  Smith  and  Clutter  (reference  7)  with  the  data  of  references  9,  11  and 
24,  coupled  to  their  own  data  for  cylindrical  elements  with  ratios  of  k/d  which 
ranged  from  0.046  to  0.4B,  concluded  that  the  critical  Rek  increased  with  k/d. 
However,  the  scatter  In  the  data  exceeded  100  percent.  In  view  of  the 
differences  among  the  various  experiments,  not  only  in  regard  to  the  type  of 
roughness  but  to  such  factors  as  the  nature  of  the  free-stream  turbulence, 
pressure  gradient,  and  as  to  how  the  critical  roughness  Reynolds  number  was 
determined,  it  was  deemed  advisable  to  assess  the  effect  of  k/d,  if  only  in  a 
limited  way,  in  the  present  Investigation  in  which  the  same  experimental 
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’'environment"  is  provided.  The  result  obtained  for  a  cylindrical  element  with  a 
k/d  of  1.0,  and  k  =  3.18  mm  at  an  xk  -  91.4  cm  is  shown  in  Figure  13-  The 

critical  Re^  has  a  value  of  450  and  is  higher  than  that  obtained  with  a 

hemispherical  element  of  comparable  k.  However,  the  values  of  325  and  450  are 
on  the  extreme  low  side  of  the  range  of  data  for  the  same  k/d  examined  by  Smith 
and  Clutter.  In  addition  to  the  effect  of  k/d  it  is  reasonable  to  expect  from  a 
stability  point  of  view  that  roughness  shape  may  play  a  role  but  the  data  in  the 
literature  is  rather  limited  and  inconclusive  on  this  point.  Klanfer  and  Owen 
(reference  25)  reanalyzed  the  data  of  reference  11  for  conical  elements  with  a 
k/d  of  0.87  and  suggested  a  critical  value  of  440  for  Re^.  Klebanoff  et  al 
(reference  9)  for  spherical  elements  In  a  lateral  row  with  various  3pacings 
obtained  an  average  value  for  the  critical  Rek  of  577,  whereas  Hall  (reference 
17),  for  isolated  spherical  elements,  obtained  critical  values  of  Rek  ranging 

from  585  to  655.  Smith  and  Clutter  (reference  7),  on  the  other  hand,  for  a 

spanwise  row  of  circular  discs  with  k/d  ranging  from  0.046  to  0.48  obtained 
values  ranging  from  100  to  500.  In  the  latter  case,  the  critical  Rek  was  that 
for  which  the  Reynolds  number  of  transition  had  decreased  to  95  percent  of  the 
value  without  roughness.  The  average  of  the  critical  Rek  they  obtained  for  when 
transition  was  considered  as  having  moved  to  the  roughness  was  600.  The  lowest 
value  of  the  critical  Rek,  for  a  three-dimensional  type  of  roughness,  was  that 
obtained  by  Hama  (reference  26)  who  used  a  spanwise  row  of  flat  equilateral 
triangular  patches,  1  mm  thick,  with  their  vertex  pointing  upstream  and  touching 
at  their  base,  and  obtained  a  value  of  45.  It  is  evident  that  there  is  a  need 
for  more  controlled  and  well  defined  experiments  in  order  to  satisfactorily 
resolve  the  trend  of  the  critical  roughness  Reynolds  number  with  roughness  size 
and  shape  for  single  roughness  elements.  An  improved  insight  into  this  behavior 
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may  well  provide  the  basis  for  an  Improved  experimental  approach  to  the  study 
of,  and  the  interpreting  of,  the  behavior  of  distributed  roughness. 

3.2  Eddy  Shedding  -  Steady  Flow 

It  has  been  known  for  some  time  (references  11,  12,  14-21)  that  associated 
with  the  behavior  of  three-dimensional  roughness  elements  there  are  periodic 
disturbances  generated  in  the  immediate  downstream  vicinity  of  the  roughness. 

The  view  stemming  primarily  from  the  flow  visualization  studies  is  that  these 
periodic  disturbances  reflect  the  presence  of  an  eddy  shedding  mechanism  that 
involves  the  generation  of  three-dimensional  vortices  having  a  "hairpin"  or 
"arch"  shape. 

The  presence  of  the  periodic  disturbances  was  readily  detected  by  hot-wire 
anemometry  in  the  present  Investigation.  They  are  illustrated  by  oscillograms 
of  the  u-fluctuatian  shown  in  figure  14.  The  oscillograms  were  obtained  down¬ 
stream  of  a  hemispherical  element  for  the  conditions,  k  -  1.7  mm,  x  -  2.54  cm, 
and  z  -  0.0  cm,  where  z  -  z-2^,  and  zk  is  the  spanwl3e  position  of  the 
roughness.  The  sweep  speed  of  5  ms  Is  noted  on  each  oscillogram.  Time 
increases  from  left  to  right,  and  a  decreasing  velocity  Is  in  a  downward 
direction.  It  is  seen  from  figure  14  that  the  frequency  increases  with  The 

distance  from  the  surface  at  which  the  oscillograms  were  obtained  was,  in 
general,  increased  with  increasing  U1 ,  reflecting  the  behavior  of  the  region  of 
periodicity  to  extend  further  from  the  surface  with  increasing  Uj  and  x.  This 
behavior  Is  consistent  with  the  flow  visualization  studies  in  which  the 
increasing  extent  of  the  "hairpin"  eddies  with  their  "heads"  moving  out  from  the 
surface,  with  increasing  and  x  was  observed. 

Wygnanski  and  Petersen  (reference  27)  showed  that  the  presence  of  vortical 
structures  and  their  behavior  inferred  from  flow  visualization  studies  may  not 
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always  be  of  dynamical  significance.  Such  structures  in  free-shear  flows  could 
be  viewed  as  a  manifestation  of  Rayleigh  instability  and  that  an  amalgamation  of 
"tagged"  particles  evidenced  by  streaklines  does  not  necessarily  coincide  with  a 
redistribution  of  vorticity.  Hama  in  an  earlier  paper  (reference  28)  had 
obtained  a  similar  result.  He  showed  that  the  rolling  up  of  a  streakline  in  a 
mixing  layer  which  had  one  fluid  at  rest  could  result  from  a  superposed  neutral 
wave,  and  therefore  could  not  necessarily  be  viewed  as  a  positive  Indication  of 
the  presence  of  vortices.  However,  references  27  and  28  do  not  necessarily 
negate  the  presence  of  vortices.  The  question  of  wave  or  vortex  may  be  viewed 
as  alternate  approaches  to  the  same  phenomenon.  Nevertheless,  it  is  desirable 
to  identify  the  physical  entity  actually  involved  since  each  view  can  provide  a 
different  mechanism  for  the  final  onset  of  turbulence. 

It  is  evident  that  the  oscillograms  of  figure  14  do  not  distinguish  between 
wave  and  vortex.  In  order  to  make  this  distinction,  the  oscillograms  shown  in 
figures  15  and  16  were  obtained  for  the  same  roughness  element  and  roughness 
position  as  for  figure  14.  In  figures  15  and  16,  as  for  figure  14,  time 
Increases  from  left  to  right,  and  a  decreasing  velocity  Is  In  a  downward  direc¬ 
tion.  However,  the  sweep  speed  was  2  ms  rather  than  5  os.  The  osoillograms  of 
figure  15  are  of  simultaneous  signals  from  two  hot-wires  separated  in  the  y 
direction  at  x  -  1.27  and  2.54  cm.  The  top  trace  In  each  oscillogram  is  the 
signal  from  the  hot-wire  that  Is  furthest  from  the  surface.  The  important 
features  are  the  phase  reversal  with  distance  from  the  surface,  which  is  clearly 
evident  with  increasing  U ^ ,  and  the  characteristic  "spike"  appearance  of  the 
fluctuation.  The  latter  permitted  the  use  of  one  hot-wire  for  detecting  the 
existence  of  a  phase  reversal  with  varying  z.  The  resulting  oscillograms 
obtained  with  U1  -  11.8  m/s,  and  at  x  -  1.27  cm.  y  -  3*18  mm  are  shown  in  figure 
16.  It  is  seen  that  phase  reversals  occur  at  +z  and  -z  with  a  spacing  which  is 


22 


AEDC -TR-87-7 


leas  than  the  diameter  of  the  roughness.  It  is  unlikely  that  a  three- 

dimensional  wave  motion  could  provide  the  appropriate  phase  reversals  and 

"spike"  appearance.  It  is  therefore  concluded,  in  accord  with  the  flow 

visualization  studies,  that  the  observed  periodicity  when  exhibiting  the 

appropriate  phase  reversals  is  indicative  of  an  eddy  shedding  frequency.  In 

figure  15,  with  U1  -  7.0  and  7.2  m/s  at  x  -  1.27  and  2.54  cm,  respectively,  it 

is  seen  that  there  is  no  discernible  phase  reversal,  although  there  is  clear 

evidence  of  a  periodic  fluctuation  at  the  inner  y  position.  At  equal  to  8.2 

and  8.3  m/s  and  x  -  1.27  and  2.54  om,  respectively,  the  phase  reversal  is  barely 

discernible  indicating  a  less  intense  vortex  than  at  the  higher  velocities.  The 

critical  velocity,  U  ,  corresponding  to  the  critical  roughness  Reynolds  number 
c 

or  325  is  8.1  m/a.  The  oscillograms  in  figure  14  also  demonstrate  the  existence 
of  periodicity  at  <  U  .  It  would  thus  appear  that  there  is  a  rapid  change 
from  wave  to  vortex  with  increasing  U1 ,  and  by  inference,  similar  behavior  with 
x  at  a  given  u  >  UQ.  The  latter  implies  an  initial  wave  motion  or  relatively 
short  wave  length  with  the  change  from  wave  to  vortex  occurring  within  a 
distance  less  than  7.5  k. 

The  dependence  of  the  frequency  on  U1  was  determined  for  the  hemispherical 
elements,  and  the  cylindrical  element  under  essentially  the  same  conditions, 
i.e.,  roughness  size  and  position,  as  for  the  determination  of  the  critical 
roughness  Reynolds  number.  The  results  obtained  with  the  hemispherical  elements 
are  shown  in  figure  17.  The  data  for  x^/k  =  359  and  393  were,  within  the 
experimental  accuracy,  considered  to  be  at  the  same  x^/k.  The  straight  lines 
are  a  least-squares  fit  to  the  data.  The  maximum  and  minimum  deviations  about 
an  average  value  of  940  Hz  shown  for  an  x^/k  ■  538  were  obtained  from  the  twelve 
oscillograms  of  figure  16.  Their  magnitude  should  not  be  regarded  as  being 
applicable  to  any  of  the  other  measurements  which  were  made  at  x  -  2.54  cm. 
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z  -  0.0  cm,  and  at  each  were  obtained  from  only  one  or  two  oscillograms  of 
the  type  shown  In  figure  1H.  The  deviations  shown,  however,  do  illustrate  the 
variability  of  the  eddy  shedding  process  and  the  lack  of  control  over  the 
frequency,  in  contrast,  for  example,  to  the  control  over  the  frequency  when 
using  a  vibrating  ribbon.  In  this  connection,  it  is  of  interest  to  note  that 
the  more  appreciable  scatter  occurs  at  frequencies  which  reflect  the  effect  of 
the  free-stream  disturbances  shown  in  figure  10,  i.e.,  1 80  Hz  and  its  harmonics. 
The  effect  of  free-stream  disturbances  in  biasing  the  maximum  amplified 
frequency  is  not  unexpected  for  a  stability  governed  phenomenon,  and  stresses 
the  importance  of  the  experimental  environment. 

The  dependence  of  the  eddy  shedding  frequency  on  U1  and  x^/k,  as  shown  in 
figure  IT,  suggests  that  the  Strouhal  behavior  can  be  nondimensionally  charac¬ 
terized  with  as  the  characteristic  velocity,  and  S*  as  the  characteristic 
length,  with  a  resulting  linear  relation  for  f,  such  that 


f 


a  C 


V 


k 


(1) 


and 


5-2C 


(2) 


The  frequency  measurements  are  shown  in  figures  18  and  19  as  suggested  by  (1) 
and  (2).  The  Strouhal  number  obtained  from  the  least-squares  fit  to  the  data  in 
figure  18,  in  accordance  with  (2),  is  0.29,  and  is  compared  in  figure  19  with 
the  actual  values  of  fC^/U^.  Implicit  in  (2)  is  the  assumption  that  Uk  varies 
linearly  with  U1  which  for  a  Blasius  flow  la  reasonably  valid  for  values  of 
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k/6*k  <  i.O.  The  range  of  k/6*k  covered  in  the  measurements  Is  shown  in  figure 
20  and  it  is  seen  that  for  Rek  greater  than  700,  the  use  of  (2)  would  tend  to 
underestimate  the  Strouhal  number.  However,  the  effect  of  k/6*k  Is  not  large, 
being  on  the  order  of  10  percent  for  the  largest  value  of  k/6#k>  Although  (2) 
appears  to  mask,  particularly  at  the  lower  Reynolds  numbers,  a  small  dependence 
on  Reynolds  number,  it  is  within  the  experimental  uncertainty,  a  reasonable 
representation  of  the  Strouhal  behavior  of  hemispherical  elements.  In  any 
event,  it  indicates  that  6*k  is  a  more  appropriate  characteristic  length  for  the 
Strouhal  number  for  single  roughness  elements  in  a  laminar  boundary  layer  than  a 
length,  such  as  k,  which  has  generally  been  used  in  characterizing  Strouhal 
behavior.  The  latter,  in  contrast,  and  as  illustrated  in  figure  21,  exhibits  a 
very  strong  dependence  on  Reynolds  number. 

The  shedding  frequencies,  and  the  resulting  Strouhal  number  obtained  with 

the  cylindrical  element,  albeit  for  only  one  size  and  position,  are  compared 

with  the  results  obtained  for  the  corresponding  hemispherical  element,  i.e.,  the 

same  k,  and  x^,  in  figures  22  and  23,  respectively.  Both  sets  of  measurements 

were  made  at  x  ■=  2.51*  cm.  It  Is  seen  that  the  shedding  frequency  for- the  cylin- 

2 

drical  element  is  also  linearly  dependent  on  U  .  However,  the  frequency  is 
lower  than  that  for  a  hemispherical  element  at  the  same  Re^  and  k/6*k.  The 
Strouhal  number  is  correspondingly  lower,  and  exhibits  a  greater  dependence  on 
Rek  at  the  lower  Reynolds  numbers,  i.e.,  for  Rek  <  800.  This  difference  In 
shedding  frequency  of  the  hemispherical  and  cylindrical  elements  Is  consistent 
with  the  greater  stability  of  the  cylindrical  roughness  in  inducing  earlier 
transition.  It  is  reasonable  to  infer  that  it  is  due  to  the  difference  in 
stability  of  the  mean  velocity  profiles  established  by  the  elements  in  their 
immediate  downstream  vicinity.  This  behavior  of  the  profiles  may  not  only  be 
reflected  in  the  greater  variation  of  ffi*k/Uk  with  Rek  for  the  cylindrical 
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element,  but  also  in  the  variation  of  the  critical  Re^  referred  to  previously  in 
Section  3. 1 . 

Data  in  the  literature  pertaining  to  the  parametric  behavior  of  the  eddy 
shedding  frequencies  is  minimal  and  limited  to  references  18  and  21.  Reference 
18  provides  data  on  the  Strouhal  behavior  for  spherical  and  cylindrical  elements 
with  the  latter  having  their  axis  in  the  spanwise  direction  rather  than  in  the 
y-directlon  as  in  the  present  study.  Reference  21,  however,  did  study  the 
Strouhal  behavior  of  hemispherical  elements.  Both  of  these  investigations,  in 
contrast  to  the  present  study,  were  conducted  in  water  channels,  and  in  neither 
investigation  are  the  results  presented  In  a  form  sufficiently  explicit  to 
permit  comparison  with  the  present  study  on  a  basis  other  than  that  shown  in 
figure  21,  Be  that  as  it  may,  the  Strouhal  behavior  for  hemispherical  elements, 
as  illustrated  in  figure  21,  is  compared  in  figure  24  with  the  Strouhal  behavior 
as  presented  in  reference  21.  The  nomenclature  of  reference  21  is  such  that  S, 
R,  xR,  and  ReR  correspond  to  fk/Uk,  k,  xk,  and  Rek,  respectively.  The  present 
measurements  are  significantly  higher  than  those  of  reference  21  over  the  range 
of  ReR.  Acarlar  and  Smith  observe  that  the  large  scatter  in  their  results  for 
Rek  >  1000  is  not  due  to  experimental  uncertainties  but  infer  from  the  occur¬ 
rence  of  multiple  peak  power  spectra  that  there  is  an  actual  instability  in  the 
shedding  behavior.  Apart  from  the  deviations  in  shedding  frequency  previously 
mentioned,  this  unstable  behavior  was  not  observed  in  the  present  study.  The 
reason  for  the  differences  between  the  present  study  and  that  of  reference  21  Is 
not  readily  apparent,  other  than  to  comment  that  in  the  present  study  the  fre¬ 
quencies  were  obtained  from  oscillograms,  whereas  in  reference  21  the  shedding 
frequency  was  determined  from  power  spectra.  It  is  perhaps  worth  commenting 
that  the  occurrence  of  multiple  peak  power  spectra  is  not  a  sufficient  condition 
for  unstable  behavior.  It  is  not  surprising  that  power  spectra  of  asymmetric 
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and  modulated  signals,  for  example,  as  shown  In  figures  14,  15  and  16  would 
exhibit  multiple  peaks  that  vary  not  only  with  but  with  y  and  x. 

The  variation  of  the  Strouhal  number,  fk/U^,  with  Rek  for  the  cylindrical 
element  with  k/d  -  1.0  is  compared  in  figure  25  with  data  faired  from  reference 
18.  The  data  selected  from  reference  18  are  for  a  spherical  element,  and  a 
cylindrical  element  with  an  aspect  ratio  of  1.0.  The  cylindrical  element  of  the 
present  study  exhibits,  within  the  experimental  uncertainty,  about  the  same 
Strouhal  numbers  for  a  comparable  range  of  Reynolds  number  as  does  the  spherical 
element.  However,  both  have  values  that  are  almost  half  of  those  for  the 
cylinder  of  reference  16.  Inasmuch  as  a  two-dimensional  roughness  element  has 
been  found  to  be  more  effective  than  a  three-dimensional  roughness  element  in 
inducing  transition,  it  is  perhaps  not  surprising  that  the  cylinder  with  its 
axis  in  the  spanwise  direction  would  exhibit  greater  instability,  if  such  is  to 
be  inferred  from  the  greater  Strouhal  number,  than  one  with  its  axis  normal  to 
the  surface.  The  effect  of  aspect  ratio,  i.e.,  length  to  diameter  ratio,  on  the 
Reynolds  number  of  transition  for  cylindrical  elements  mounted  on  the  surface 
with  their  axis  in  the  spanwise  direction  was  examined  by  Norman  {reference  19). 
The  Reynolds  number  of  transition,  a3  observed  at  a  fixed  distance  downstream 
from  the  element,  was  found  to  decrease  with  increasing  aspect  ratio,  over  a 
range  of  aspect  ratios  from  1  to  8,  and  at  an  aspect  ratio  of  8  the  Reynolds 
number  of  transition  was  actually  lower  than  for  a  two-dimensional  element. 
Furuya  and  Miyata  (reference  18),  on  the  other  hand,  Investigated  the  effect  of 
aspect  ratio  on  Strouhal  numbers  for  cylindrical  elements  with  their  axis  in  the 
spanwise  direction.  They  found,  over  a  range  of  aspect  ratios  from  1  to  6  that 
the  Strouhal  numbers  decreased  with  increasing  aspect  ratio.  The  f orementloned 
results  from  references  18  and  19  lead  to  the  inference  of  decreasing  Strouhal 
numbers  with  decreasing  critical  roughness  Reynolds  numbers  for  cylindrical 
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elements  with  their  axis  in  the  spanwise  direction.  This  trend  of  Strouhal 
number  with  critical  roughness  Reynolds  number  is  opposite  in  direction  from 
that  observed  in  the  comparison  of  the  behavior  of  hemispherical  elements,  and  a 
cylindrical  element  in  the  present  study.  The  issue  is  complicated  by  the 
evidence  (reference  10)  that  two-dimensional  roughness  elements  in  inducing 
earlier  transition  do  not  exhibit  the  critical  behavior  involving  the  shedding 
of  discrete  vortices.  Roughness  shape,  as  well  as  aspect  ratio,  appear  to  play 
a  role,  and  It  may  well  be  that  the  Strouhal  number  is  not  a  sufficient  criter¬ 
ion  for  Indicating  the  degree  of  instability.  It  is  apparent  that  these  aspects 
warrant  further  3tudy.  Their  resolution  Is  amenable  to  experiment,  and  would  go 
a  long  way  toward  providing  a  better  understanding  of  the  mechanism  by  which 
discrete  roughness  elements  induce  earlier  transition  to  turbulent  flow. 

3.3  Eddy  Shedding  -  Unsteady  Flow 
One  would  readily  infer  from  the  discussion  in  Section  3.2  that  the 
"hairpin"  eddies  are  intrinsic  to  the  transition  process  induced  by  a  three- 
dimensional  roughness  element.  However,  the  view  existing  in  the  literature 
(references  1*1,  19,  20,  29,  30)  is  that  such  eddies  are  not  related,  or  at  be3t 
are  only  indirectly  related,  to  the  transition  process.  It  was  therefore  deemed 
desirable  to  attempt  to  clarify  the  role  they  do  play.  To  this  end,  the  effect 
of  an  oscillatory  free-stream  was  Investigated.  The  free-stream  velocity,  , 
was  oscillated  sinusoidally  with  a  small  amplitude,  AU^ ,  about  some  mean  value, 

U  ,  and  at  a  sufficiently  low  frequency,  n,  for  a  quasi-steady  assumption  to  be 
valid.  That  is. 


Uj (t)  -  U1  (1  +  A  sin  ut) 


(3) 
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where  A  =  AU^/l^,  and  u  =  2irn.  Of  particular  Interest  is  the  effect  of  the 

variation  of  U,  relative  to  U  . 

1  c 

The  response  of  the  laminar  boundary  layer  without  roughness  to  the  free- 
stream.  oscillation  is  illustrated  in  figure  26.  Figure  26  shows  the 
distributions  across  the  boundary  layer  of  the  oscillating  velocity  component, 
AU,  and  the  phase  of  the  oscillation,  $.  A  positive  angle  represents  an  advance 
in  phase  relative  to  the  free-stream  oscillation.  The  measurements  were  made  at 
x  =  94  cm  with  ^  -  6.8  m/s  and  n  =*  1.0  Hz.  The  resulting  frequency  parameter, 
x^,  given  by  xu/i^ ,  was  0.87,  and  A  was  0.12.  The  trends  exhibited  by  the 
distributions  are  similar  to  previous  results  (reference  31).  The  mean  velocity 
distribution  for  the  unsteady  boundary  layer  without  roughness  is  compared  with 
its  steady  flow  counterpart  in  figure  27.  Also  shown  in  figure  27  is  a  com¬ 
parison  of  the  mean  velocity  profiles  for  steady  and  unsteady  flow  downstream  of 
a  hemispherical  roughness  element.  In  the  latter  case,  the  measurements  were 
made  on  the  centerline  at  2. 54  cm  downstream  of  a  hemispherical  element  which 
had  a  nominal  diameter  of  3.18  mm,  an  effective  height  of  1.7  mm  and  was  posi¬ 
tioned  at  xk  =  91.4  cm.  The  free-stream  velocities  were  slightly  above  critical 
being  8.3  and  8.2  m/s  for  the  steady  and  unsteady  flow  respectively.  The 

oscillating  frequency  wa3  1  Hz  with  A  -  0.13.  It  can  be  concluded  from  the 

agreement  of  the  disturbances  for  steady  and  unsteady  flow  that  the  oscillation 
behaves  linearly.  In  addition,  It  is  of  interest  to  note  that  the  mean  velocity 

distribution  at  x  -  14.9  k,  at  about  the  critical  speed,  la  insensitive  to  the 

time  dependence  of  the  eddy  generating  process  associated  with  the  unsteady 
flow. 

The  effect  of  an  oscillating  free-stream  on  the  eddy  shedding  process  was 
studied  with  the  hemispherical  element  positioned  at  xk  ■*  91.4  cm  with  a  k  -  1.7 
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ram.  Simultaneous  signals  of  u  in  the  boundary  layer  and  of  in  the  free- 
stream  with  free-stream  oscillation  frequencies  of  1  and  2  Hz,  at  various  and 
associated  AU1 ,  were  monitored  by  strip-chart  recording  at  z  -  0.0  cm  and  x  ■» 

2. 54,  30.5,  and  71.0  cm.  Figures  26-37  show  the  recordings  that  were  obtained. 
The  upper  trace  in  each  recording  is  the  free-atream  oscillating  velocity,  and 
the  lower  trace  Is  the  u-fluctuation  in  the  boundary  layer.  Time  increases  from 
left  to  right,  and  a  decreasing  velocity  is  in  a  downward  direction.  Figures  28 
and  29  show  the  recordings  at  x  -  2.5*1  cm  for  oscillation  frequencies  of  1  and  2 
Hz  respectively.  The  position  in  the  boundary  layer  was  at  y  «*  3-18  mm.  The 
eddy  shedding  is  characterized  by  intense  fluctuations  In  the  direction  of  lower 
velocity  superposed  on  the  primary  oscillation.  This  type  of  signal  la 
consistent  with  the  passage  of  "hairpin"  eddies  with  "legs"  that  are  rotating  as 
indicated  In  Section  3.2.  It  should  be  noted  that  the  asymmetry  of  eddy 
shedding  relative  to  the  free-stream  oscillation  as  shown  in  figures  28  and  29 
is  not  a  real  effect.  It  should  not  be  interpreted,  for  example,  as  arising 
from  a  difference  in  the  stability  of  an  accelerating  flow  via  a  vis  a 
decelerating  flow.  It  arises  from  an  inadvertent  low-frequency  cut-off  of  2  Hz 
in  the  hot-wire  electronics  for  both  probes  which  Introduced  a  60  degree  and  a 
40  degree  phase  advance  relative  to  the  phase  of  the  eddies  for  oscillation 
frequencies  of  1  and  2  Hz  respectively.  In  addition,  the  low-amplitude  ripple 
of  20  Hz  in  the  boundary  layer  signal  is  a  consequence  of  the  vibration  induced 
by  the  fan  motor  which  was  previously  referred  to  in  Section  2.3  and  which 
manifests  itself  because  of  the  velocity  gradient  in  the  boundary  layer.  In  any 
event,  these  extraneous  effects  do  not  mitigate  the  essential  features  of  the 
effect  of  unsteady  flow  on  the  behavior  of  the  roughness. 

It  is  seen  from  figure  28  that  for  a  steady  flow  at  U1  -  8.5  m/s  (Rek  - 
337)  that  eddies  are  shed  continuously.  When  the  flow  is  oscillated  at  a  value 
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of  U1  which  coupled  to  the  associated  A^  is  below  some  value  of  U1  then  no 

eddies  are  observed  as  in  the  recording  for  ■»  5.8  m/s  with  AU.^  »  0.64  m/s. 

However,  eddies,  although  of  relatively  lower  intensity,  are  observed  over  part 
of  the  cycle  in  the  recording  for  ■  6.9  m/s  and  AU^  =  0.84  m/s.  As  is 

further  increased,  the  duration,  AT,  of  the  eddies  increases,  and  at  U1  -  9.7 

m/s  with  AU1  -  1.38  m/s  eddies  are  shed  continuously  over  the  cycle.  A 
recording  with  a  much  faster  chart  speed,  for  the  latter  condition,  is  also 
shown  in  order  to  more  clearly  illustrate  the  nature  of  the  fluctuations  in  the 
recording  and  their  similarity  to  the  fluctuations  for  the  steady  flow 
condition.  If  an  appropriate  allowance  is  made  for  the  aforementioned  phase 
advance,  the  variation  of  the  amplitude  and  frequency  of  the  fluctuation  with 
is  consistent  with  quasi-steady  behavior.  The  recordings  in  figure  29  3how 
results  similar  to  figure  28.  The  smaller  amplitude  at  the  various  U1  with 
n  -  2  Hz,  as  compared  to  n  -  1  Hz,  is  a  consequence  of  the  time  constant  of  the 
Unsteady  Wind  Tunnel.  However,  it  is  evident  that  except  for  the  effect  the 
amplitude  may  have  on  AT  the  behavior  scales  with  the  frequency  n. 

It  appears  from  the  recordings  In  figures  28  and  29  that  the  onset  of  the 
"hairpin"  eddies  is  associated  with  a  critical  value,  U1c,  of  the  oscillatory 
free-stream.  It  can  be  readily  shown  from  equation  (3)  that 

AT  -  ^  [w  -  2  sin"1  (  ^  '■■■>]  (1° 

and 

uic  “  AU1  3in  [  [l  '  AT)  P  +  U1  (5> 
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It  is  seen  from  equation  (4)  that  the  scaling  of  AT  with  frequency  n  is  to  be 
expected  and  that  the  effect  of  amplitude  is  to  change  the  slope  of  the 
variations  of  AT  with  U1  about  AT  »  0.5.  It  is  of  interest  to  compare  the 
values  of  U1c  with  the  critical  velocity,  Uq,  corresponding  to  the  critical 
roughness  Reynolds  number  of  transition  observed  In  Section  3.1.  The  values  of 
UlQ  obtained  from  equation  (5)  for  the  various  recordings  In  figures  28  and  29, 
using  the  measured  values  of  AT,  averaged  over  a  number  of  cycles,  are 
summarized  In  Table  I. 

Table  I 


U1 (m/s) 

AU1 (m/s) 

n(Hz) 

AT(s) 

Ule(m/s) 

6.9 

0.84 

1.0 

0.19 

7.6 

7.6 

0.98 

1.0 

0.36 

8.0 

8.7 

1.20 

1.0 

0.68 

8.1 

7.1 

0.58 

2.0 

0.09 

7.6 

7.6 

0.64 

2.0 

0.19 

7.9 

8.7 

0,78 

2.0 

0.39 

8.1 

No  attempt  was  made  to  establish  with  specificity  the  conditions  for  AT  <*  1 ,  and 

consequently  the  associated  value  of  U.  is  not  included  In  the  table.  The 

1  c 

consideration  of  a  critical  value  of  velocity,  U^c,  from  the  recordings  in 

figures  28  and  29  is  not  Inconsistent  with  the  discussion  in  Section  3.2,  i.e., 

the  presence  of  periodic  disturbances  at  values  of  U,  <  U  and  that  there  is  a 

l  c 

continuous  transition  from  wave  to  vortex  with  increasing  Uj.  It  is  not 

inconsistent  because  of  the  inference  from  figures  28  and  29  that  disturbances 

which  may  exist  for  U.  <  U,  ,  are  restricted  in  their  extent  from  the  surface  to 

l  1  c 

y  <  3-18  ram.  Concomitantly,  at  a  given  and  AU^ ,  Ulc  is  not  a  uniquely 

critical  velocity.  It  may  vary  with  x  and  y,  depending  on  the  topology  of  the 
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vortical  structures.  Nevertheless,  the  correspondence  of  the  values  of  and 

U  «  8.3  m/s  (Re,  -  325)  is  consistent  with  the  conclusion  that  the  eddy 
c  K 

generation  is  a  necessary  condition  for  the  Induced  transition.  The  difference 

in  the  value  of  U  -  8.3  m/s  from  the  previously  stated  value  of  8.1  m/s  used  in 
c 

connection  with  the  discussion  of  figures  14  and  i5  arises  from  the  variation  in 

v. 

The  question  as  to  the  intrinsic  importance  of  the  eddies  to  the  transition 
process  and  the  turbulence  structure  is  addressed  further  by  the  strip-chart 
recordings  shown  in  figures  30-37.  These  were  obtained  over  a  range  of  and 
AU1  comparable  to  those  in  figures  28  and  29.  The  various  positions  in  the 
boundary  layer  and  frequencies  of  the  oscillating  free  stream  velocity  for  which 
the  recordings  were  obtained  are  summarized  in  Table  II. 


Table  II 

Figure  No.  x(cm)  y(mm)  n(Hz) 

30  30.5  0.61  1.0 

31  30.5  0.61  2.0 

32  30.5  6.35  1.0 

33  30.5  6.35  2.0 

34  61.0  0.61  1.0 

35  61.0  0.61  2.0 

36  61.0  6.35  1.0 

37  61.0  6.35  2.0 

As  for  figures  28  and  29,  time  Increases  from  left  to  right,  and  a  decreasing 
velocity  is  in  a  downward  direction.  The  recordings  reveal  a  number  of 
interesting  features.  It  is  seen  that  the  eddies  which  are  clearly  observable 
at  x  «  2.54  cm  have,  by  x  -  30.5  cm,  induced  or  evolved  into  a  turbulent  spot. 


33 


AEDC-TR-S7-7 


In  this  connection,  it  should  be  noted  that  the  Inadvertent  low  frequency  cut¬ 
off  of  2  Hz,  referred  to  previously,  has  distorted  the  signal  associated  with 
the  "spot"  in  the  recordings  with  an  oscillating  free-stream.  More 
specifically,  It  is  responsible  for  the  decreasing  velocity  associated  with  the 
spot  passage,  and  the  minimum  and  rapid  rise  following  the  passage  of  the  spot 
as  seen,  for  example,  in  the  recording  for  ■»  8.2  m/s  in  figure  30.  Apart 
from  this  extraneous  effect  the  "spots"  as  shown  in  the  various  recordings  of 
figures  30-37  exhibit  the  characteristic  features  generally  associated  with  the 
passage  of  a  turbulent  "spot"  (references  32-34).  These  are  an  abrupt  step-up 
in  velocity  near  the  surface,  and  abrupt  step-down  in  velocity  in  the  outer 
region  of  the  boundary  layer,  and  an  increase  in  "spot"  duration  with  x 
resulting  from  the  difference  in  velocity  at  the  leading  and  trailing  interfaces 
of  the  "spot".  The  changes  in  velocity  associated  with  the  "spot"  passage 
reflect  the  difference  between  the  turbulent  velocity  profile  within  the  "spot", 
and  the  laminar  velocity  profile  in  the  flow  surrounding  the  "spot".  However,  a 
cautionary  comment  should  be  introduced  in  regard  to  the  "spot"  shape.  It  is 
not  likely  that  a  turbulent  "spot"  that  originates  from  a  finite  train  of 
eddies,  as  in  the  present  case,  will  be  similar  in  shape  to  the  "spot"  which 
originates  from  a  momentary  disturbance  as  in  references  32-34. 

For  a  steady  velocity,  which  is  near  the  critical  speed,  Uc,  random 
turbulent  "spots"  form  downstream  of  the  roughness  element.  This  is  evident 
from  each  of  the  recordings  in  figures  30-37  for  which  -  8.2  or  8.3  ®/fli  and 
n  »  0.0  Hz.  The  corresponding  values  of  Rek  are  31 2  and  319  for  U1  -  8.2  and 
8.3  m/s  respectively.  It  appears  that  for  the  steady  flow  condition  at  x  -  61.0 
cm  that  the  "spots"  are  not  only  much  more  developed  in  their  turbulence 
"structure"  as  compared  to  x  -  30.5  cm  but  may  have  Increased  In  number.  It  is 
not  unexpected  that  the  "spots"  will  grow  in  size  as  they  proceed  downstream, 
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however  it  is  not  clear  whether  the  aforementioned  differences  between  the 
"spots"  at  it  »  61.0  cm,  and  B0.5  cm  are  real  or  apparent.  Thi3  aspect  was  not 
investigated  further,  however,  the  ’"spot"  formation  is  extremely  sensitive  to 
U1 ,  and  it  may  be  that  the  differences  reflect  the  imprecision  In  U^. 

It  is  seen  from  the  steady  flow  condition  in  figures  28  and  29,  at  U1  »  8.5 
m/s  (Rek  =  337)  that  the  Intensity  of  the  eddies  is  highly  modulated.  This 
behavior  provides  an  explanation  for  the  random  occurrence  of  turbulent  "spots" 
downstream  of  the  roughness  element,  although  there  is  a  continuous  generation 
of  eddies  by  the  roughness  element.  It  Is  that  the  eddy  shedding  process,  not 
being  well  controlled,  and  being  highly  modulated,  eddies  of  relatively  weak 
intensity  may  damp,  and  not  all  eddies  participate  In  the  evolution  of  the 
turbulent  "spot".  This  point  of  view  Is  supported  by  the  behavior  with  an 
oscillating  free-stream  at  the  lower  values  of  U^m),  where  U^m)  =  VJ1  +  &U1  . 

In  the  recordings  for  (m)  =>  7.86,  7.74,  and  7-68  m/s  in  figures  30,  34  and  35, 
respectively,  no  turbulent  spots  are  observed,  although  eddy  shedding  is 
observed  in  the  recordings  for  U^m)  «■  7.74  and  7-68  m/s  in  figures  28  and  29. 

In  addition  the  recordings  at  U^m)  -  8.22,  8.10,  and  7.98  m/s  in  figures  30, 

34,  and  36,  respectively,  show  evidence  of  intermittency,  i.e.  "spots"  do  not 
necessarily  occur  every  cycle.  It  should  be  noted  that  the  above  mentioned 
conclusion  that  the  recordings  in  figure  28  and  29  at  (m)  •  7.74  and  7*68  m/s 
Illustrate  eddy  shedding  is  inferred  from  the  decreasing  velocity  fluctuation 
associated  with  the  disturbance,  and  that  no  attempt  was  made  to  specifically 
Identify  the  vortical  nature  of  the  disturbance.  In  this  connection.  Hall 
(reference  17)  in  his  smoke  visualization  study  of  the  flow  behind  a  sphere 
observed  that  although  strong  periodic  laminar  vortices  were  formed  at 
subcritical  Reynolds  numbers,  no  turbulence  occurred  downstream,  and  the 
vortices  decayed  as  they  propagated  downstream.  Although  Norman  did  not 
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explicitly  address  this  question  in  reference  19,  Morkovin  (reference  35)  has 
commented  that  Norman  did  observe  similar  behavior  behind  a  fence-type  roughness 
element  with  an  intensity  level  as  high  as  4  percent. 

The  presence  of  random  turbulent  "spots"  downstream  of  the  roughness,  at 
about  the  critical  velocity,  provides  a  different  view  than  the  flow 
visualization  studies  in  regard  to  the  detailed  behavior  of  the  transitional 
region  where  transition  moves  rapidly  upstream  as  3hown  in  figure  12.  The 
china-clay  studies  of  references  11  and  16,  and  the  long-exposure  smoke 
photographs  of  reference  14,  have  led  to  the  view  that  the  rapid  upstream 
movement  of  transition  is  due  to  the  vertex  of  a  turbulent  wedge  that  moves 
rapidly  forward  with  a  slight  increase  in  roughness  Reynolds  number.  The  view 
from  the  present  study  is  that  this  transitional  region  consists  of  turbulent 
"spots"  which  increase  in  number  with  a  slight  increase  in  roughness  Reynolds 
number.  As  they  move  downstream,  the  "spots"  grow  laterally,  overtake  one 
another,  and  sweep  out  a  wedge-like  boundary  which  gives  the  impression  of 
forward-moving  turbulent  wedge  In  flow  visualization  studies  with  a  relatively 
long-time  exposure.  However,  at  a  value  of  Rek,  sufficiently  above  the  critical 
roughness  Reynolds  number,  the  "spots"  may  have  so  Increased  in  number,  and/or 
the  eddies  being  3hed  may  have  so  Increased  in  intensity,  that  a  turbulent  wedge 
condition  consistent  with  the  flow  visualization  studies  has  been  established. 

It  is  shown  in  figures  28-37  that  as  U1  increases,  the  duration  of  the 
turbulent  "spots”  as  well  as  the  duration  of  eddy  shedding  Increases.  It  was  of 
particular  interest  to  examine  how  the  duration  of  "spots"  and  eddies  compared. 
For  this  purpose,  the  variation  of  nAT  with  at  x  ■  2.54,  30.5,  and  61.0  cm  is 
shown  in  figure  38.  Within  the  experimental  accuracy,  the  small  differences 
between  the  duration  of  eddy  shedding  at  1  Hz  and  2  Hz,  at  x  ■»  2.54  cm,  may  be 
real  and  reflect  the  effect  of  amplitude  as  given  by  equation  (4).  This  effect 


36 


AEDC-TR-87-7 


may  also  influence  the  scatter  in  the  data  at  x  -  30.5  era.  However,  faired 
curves  are  drawn  through  the  data  in  order  to  facilitate  the  comparison.  It  Is 
seen  that  the  degree  of  correlation  increases  with  U1  .  For  example,  at  >■  8 
m/s,  the  duration  of  the  turbulent  "spot"  at  x  -  30.5  cm  is  0.73  the  duration  of 
eddy  shedding.  At  =9  m/3,  the  duration  of  the  turbulent  "spot"  at  x  -  30.5 
cm  is  0.9  the  duration  of  eddy  shedding.  This  behavior  supports  the  view 
previously  expressed,  i.e.,  that  with  increasing  Re^  the  intensity  or  the  eddies 
may  Increase,  and  a  greater  number  nay  participate  in  the  development  of  the 
"spot".  At  x  =  61.0  era,  the  duration  of  the  "spot"  at  -  8.0  and  9.0  m/s  is 

greater  than  that  at  x  -  30.5  cm.  This  is  not  unexpected  since  the  velocity  of 
the  leading  edge  of  the  "spot"  is  much  larger  than  the  velocity  of  the  trailing 
edge,  and  a  small  difference  in  AT  can  represent  a  large  difference  in  the  size 
of  the  "spot".  The  high  degree  of  correlation  of  the  duration  of  turbulent 
"spots"  with  the  duration  of  eddy  shedding,  coupled  to  the  fact  that  when  the 
eddy  shedding  is  periodic  with  the  frequency  of  the  free-stream  oscillation  the 
turbulent  "spots"  also  became  periodic,  warrant  the  conclusion  that  not  only  is 
the  eddy  shedding  essential  to  the  transition  process,  but  the  eddieq,  although 
modified  in  form,  are  intrinsic  to  the  vortical  structure  of  the  turbulence. 

3.4  Instability  and  Induced  Transition 
In  their  study  of  the  transition  process  for  a  Blaslus  flow  on  a  flat 
plate,  Klebanoff  et  al.  (reference  36)  drew  an  analogy  between  the  secondary 
inflectional  instability  occurring  In  the  later  stage  of  transition,  and  the 
behavior  of  a  three-dimensional  roughness  element.  The  "hairpin"  eddies 
described  in  Sections  3.2  and  3-3  are  markedly  similar  to  the  eddies  described 
in  reference  36  in  connection  with  the  secondary  instability.  However,  the  eddy 
shedding  frequency  associated  with  the  roughness,  as  shown  in  figure  17,  varies 
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2  3/2 

with  U1  rather  than  with  as  in  reference  36,  and  the  analogy  may  best  be 

regarded  as  qualitative. 

The  nature  of  the  flow  about  a  three-dimensional  roughness  element  which 
results  in  an  Inflectional  velocity  profile  in  the  Immediate  downstream  vicinity 
of  the  roughness  has  been  well  documented  (references  11-21).  The  essential 
features  disclosed  primarily  by  flow  visualization  studies  are  that  two  sets  of 
vortices  are  established  at  a  subcritlcal  Reynolds  number.  One  is  a  closely 
spaced  pair  of  spiral  filaments  which  form  in  the  near  wake,  3piral  upward  at 
the  rear  of  the  roughness,  and  at  the  height  of  the  roughness  turn  and  trail 
downstream.  The  other  is  a  "horseshoe"-shaped  vortex,  close  to  the  surface, 
which  wraps  around  the  front  of  the  roughness  forming  a  pair  of  streamwlse 
vortices  which  also  extend  downstream.  This  array  of  vortices  which  results  In 
an  Inflectional  profile  at  z  ■  0.0  cm  has  been  observed  for  a  variety  of  three- 
dimensional  roughness  geometries  such  as  spherical,  cylindrical,  hemispherical, 
conical,  and  rectangular.  The  flow  visualization  studies  have  also  yielded  the 
observation  that  as  the  free-stream  velocity  is  increased,  a  waviness  appears  on 
the  vortex  filaments,  and  this  waviness  increases  with  increasing  velocity  until 
the  filaments  "roll-up"  into  "hairpin"  eddies. 

The  qualitative  behavior  revealed  by  the  flow  visualization  studies 
manifests  itself  in  the  spanwi3e  distributions  of  U  and  u‘  as  shown  in  figure  39 
and  MO  respectively.  The  measurements  were  made  at  a  distance  of  1M.9  roughness 
heights  downstream  of  the  roughness  elements  and  for  three  different  unit 
Reynolds  numbers.  The  behavior  of  the  frequency  of  eddy  shedding  discussed  in 
Section  3.2  suggested  that  it  would  be  more  appropriate  to  carry  out  the 
boundary  layer  surveys  at  a  given  U^/v  rather  than  .  Consequently,  was 
adjusted  to  compensate  for  changes  in  v,  and  the  value  of  noted  in  figure  39 
is  the  average  value  for  the  various  distributions.  The  unit  Reynolds  numbers 


38 


AE0C-TR-B7-7 


selected  were  such  that  one  was  below  the  critical  roughness  Reynolds  number, 
one  was  at  or  slightly  above  critical,  and  one  was  significantly  above  critical. 
In  this  connection,  it  should  be  noted  that  the  value  of  U^/v  corresponding  to  a 
critical  roughness  Reynolds  number  of  325  is  *1.76  x  10^  (m  ' ) .  At  the 
subcritical  Reynolds  number  the  spanwise  distributions  of  U  clearly  demonstrate 
the  presence  of  the  two  pairs  of  streamwise  vortices  observed  in  the  flow 
visualization  studies.  The  pair  of  vortex  filaments  rotate  in  such  a  direction 
as  to  transfer  momentum  away  from  the  plate  at  z  =  0.0  cm,  and  the  ’’horseshoe1' 
vortex  rotates  so  as  to  transfer  momentum  towards  the  plate  at  z  »  +  0.5d.  The 
spanwise  distribution  of  u'  at  the  3ubcritical  Reynolds  number  apparently 
reflects  a  low  level  unsteadiness  in  the  two  sets  of  stationary  vortices.  At 
the  higher  Reynolds  numbers  the  peak  in  the  spanwise  distributions  of  u'  in 
figure  40  closely  indicates  the  presence  of  the  "hairpin"  vortices.  With 
increasing  Reynolds  number  the  peak  increases  in  intensity  and  extends  further 
from  the  surface  consistent  with  the  behavior  of  the  "hairpin"  eddy  as  discussed 
in  Sections  3.2  and  3.3.  It  can  also  be  inferred,  within  the  accuracy  of  the 
spanwise  distributions,  that  the  spanwise  extent  of  the  eddy  is  on  the  order  of 
z/d  >  +  0.5,  centered  about  z  -  0.0  cm,  and  that  the  disturbed  boundary  layer 
can  be  separated  into  two  regions.  An  outer  region,  y/k  ^1.0,  where  the 
"hairpin"  eddy  is  dominant,  and  an  inner  region,  y/k  <  1.0,  where  the  two  sets 
of  vortices  alluded  to  at  the  subcritical  Reynolds  number,  still  play  a  role. 

The  latter  conclusion  is  supported  by  the  double  peaks  in  the  spanwise 
distribution  of  u1 ,  close  to  the  surface  at  the  Reynolds  numbers  above  critical, 
which  have  the  same  spanwise  spacing  a8  the  peaks  at  the  subcritical  Reynolds 
number.  The  intensity  of  these  peaks  increases  markedly  from  about  0.15  percent 
at  the  subcritical  Reynolds  number  to  about  6  percent  at  the  highest  Reynolds 
number.  However,  this  increase  in  intensity  is  not  necessarily  inconsistent 
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with  the  presence  of  the  two  sets  of  vortices.  The  increase  in  Intensity  may 
reflect  an  increase  in  concentration  of  the  vortices  with  Increasing  velocity  as 
well  as  a  degree  of  unsteadiness  resulting  from  the  interaction  with  the 
"hairpin"  vortices. 

The  spanwise  distributions  of  U  at  the  higher  Reynolds  number,  in  the 
presence  of  eddy  shedding,  are  similar  in  appearance  to  the  distributions  at  the 
subcritlcal  Reynolds  number  except  that  they  exhibit  a  greater  spanwise 
variation  in  U  which  extends  over  a  greater  region  of  the  boundary  layer.  Be 
that  as  it  may,  the  spanwise  distributions  of  U  can  also  be  interpreted  as  being 
consistent  with  the  separation  of  the  disturbed  boundary  layer  into  outer  and 
inner  regions  as  outlined  above.  However,  the  "hairpin"  eddies  do  have  the  same 
sense  of  rotation  as  the  set  of  spiral  vortices,  and  their  orientation  may  be 
such  as  to  make  the  above  interpretation  somewhat  tenuous.  Nevertheless,  it  is 
surprising  that  the  change  from  a  stationary  set  of  spiral  vortices  to  a  tine 
dependent  "hairpin"  eddy  shedding,  deduced  from  the  flow  visualization  studies, 
does  not  have  a  greater  effect  on  the  spanwise  distributions  of  U  and  u’  near 
the  surface. 

In  order  to  examine  the  instability  question  further,  as  well  a3  to  examine 
the  transition  to  a  fully  developed  turbulent  boundary,  profiles  of  U  and  u' , 
across  the  boundary  layer,  were  measured  on  the  centerline  for  different 
downstream  positions  ranging  from  x  -  1.27  to  30.5  cm.  The  resulting 
distributions  for  U  and  u*  are  shown  in  figures  41(a)-(g)  and  42(a)-(g)  1 

respectively.  The  distributions  at  each  x-positlon,  except  those  at  x  ■  30.5  cm 
were  measured  at  five  different  unit  Reynolds  numbers,  three  of  which  were  the 
same  as  for  the  spanwise  distributions  of  figures  39  and  *10,  one  was  an 
additional  unit  Reynolds  number  below  critical,  and  one  an  additional  unit 
Reynolds  number  above  critical. 
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The  mean  velocity  profiles  at  x  =  1.27  cm  in  figure  41(a)  are  sufficiently 

inflectional  to  sustain  an  incipient  wave  type  of  instability  which  develops 

into  a  rolled-up  vortex,  and  eddy  generation  at  and  above  the  critical  roughness 

Reynolds  number.  The  Inflectional  nature  of  the  mean  velocity  profiles  in  the 

immediate  downstream  vicinity  of  the  roughness  at  x  =  1.27  cm  (7.5k),  and 

z  -  0.0  cm,  13  illustrated  in  figures  43(a)  and  43(b)  in  which  they  are  compared 

with  inflectional  profiles  for  which  stability  calculations  have  been  carried 

out.  In  figure  43(a)  the  mean  velocity  profiles  at  unit  Reynolds  up  to  and 

Including  the  unit  Reynolds  number  which  is  slightly  above  critical  are  compared 

with  the  Lin  profile  (reference  37)  for  a  two-dimensional  mixing  layer  in  which 

one  fluid  is  at  rest.  Similarly,  in  figure  43(b)  the  same  profiles  are  compared 

with  the  Hartree  profile  (reference  38)  for  a  two-dimensional  separated  layer. 

The  measurements  of  mean  velocity,  close  to  the  surface,  and  which  are  nearly 

constant,  reflect  the  behavior  of  a  hot-wire  in  a  separated  region,  and  the 

momentum  thickness,  8,  was  obtained  by  extrapolating  the  distributions  smoothly 

to  zero  at  the  surface.  The  boundary  layer  profiles  were  matched  to  the  mixing 

layer  distribution  at  U/U^  -  0.5,  and  it  is  seen  from  the  comparison,,  as  well  as 

from  the  comparison  in  figure  43(b),  that  the  inflection  is  sufficiently  close 

to  the  wall  that  the  instability  may  be  regarded  as  being  more  characteristic  of 

a  wall-bounded  flow  rather  than  a  free-ahear  layer,  but  one  which  is  specific  to 

the  geometry  of  the  roughness,  and  the  Reynolds  number.  This  is  supported  by 

2 

the  variation  of  frequency  with  U T  ,  at  least  for  the  size  of  roughness  elements 

herein  investigated,  whereas  a  free-shear  layer  would  imply  a  variation  with 
3/2 

.  In  addition,  the  mixing  layer,  and  the  separated  laminar  boundary  layer 
have  their  inflection  at  U/U^  -  0.5.  The  presence  of  the  separated  region  tends 
to  obscure  the  inflection  point  in  the  boundary  layer  profiles.  However,  the 
corresponding  distributions  of  u*  at  x  -  1.27  cm  in  figure  42(a)  indicate  a 
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maximum  in  u*  at  y/k  =  1.34.  It  is  reasonable  to  assume  that  this  position,  for 
a  resulting  instability,  is  the  point  of  inflection.  The  corresponding  value  of 
U/U1 ,  averaged  over  the  three  unit  Reynolds  numbers,  is  0.34  which  at 
U1  •  8.3  m/a  yields  6.5  mm  for  the  wavelength  of  the  instability.  From  this 
point  of  view,  the  instability  at  the  critical  Reynolds  number  i3  strongly 
nonlinear  with  u'/U^  of  about  6  percent,  and  the  "hairpin”  eddy  developing 
within  one  or  two  wavelengths.  This  rapid  development  of  the  "hairpin"  eddy  is 
not  inconsistent  with  an  inflectional  instability  when  it  is  compared  to  the  two 
to  three  wavelengths,  observed  In  references  36  and  39,  over  which  the  nonlinear 
behavior  of  three-dimensional  Tollmien-Schlichting  waves  resulted  In  the 
generation  of  "hairpin"  eddies.  However,  in  the  latter  case  the  wavelengths 
were  very  much  longer. 

The  distributions  of  u'  indicate  that  at  the  aubcritical  Reynolds  numbers 

the  observed  instability  may  behave  linearly,  i.e.  governed  by  linear  stability 

5  -1 

theory.  At  U1  /v  -  4.13  x  10  (m  )  the  maximum  In  u1/^  increases  from  t.O 
percent  at  x  =  1.27  cm  to  1.8  percent  at  x  -  2.54  cm,  and  then  damps  with 

5  -1 

increasing  distance  downstream.  At  the  lower  Reynolds  number  of  3.54  x  10  Cm  ) 
the  maximum  in  u'/U1  at  x  =  1.27  cm  is  0.65  percent,  and  the  nature  of  the 
distribution  is  such  as  to  Indicate  that  this  represents  a  small  degree  of 
amplified  disturbance  above  the  background  disturbances.  However,  in  contrast 
to  the  higher  Reynolds  number,  it  exhibits  no  further  growth  and  damps  with 
distances  further  downstream.  It  was  shown  in  reference  39  that  a  relatively 
weak  three-dimensional  Tollmien-Schlichting  wave  with  an  intensity  level  of  2 
percent  may  still  behave  linearly.  It  is  therefore  reasonable  to  infer  that  the 
Instabilities  at  the  subcritlcal  Reynolds  numbers  are  also  behaving  linearly. 

In  this  connection,  it  is  of  interest,  despite  the  three-dimensionality  of  the 
instability,  to  compare  the  unstable  frequency  observed  at 
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5  “1 

U^v  -  M.13  x  10  (m  )  with  the  two-dimensional  stability  calculations  of 

references  *40  and  41  for  the  separated  boundary  layer  and  mixing  layer 

respectively.  The  pertinent  nondimens lonal  frequency  parameters,  2¥fv/U  2  and 

2vf0/U1  are  68  x  10  ^  for  the  separated  layer,  and  0.20  for  the  mixing  layer. 

The  corresponding  Reynolds  numbers  based  on  displacement  thickness  and  momentum 

thickness,  Rg#  and  R^,  are  1170  and  300.  At  2Trfv/U  1 2  -  68  x  10_5,  Branch  II  of 

the  stability  diagram  for  the  separated  layer  is  at  R  =  1,000.  At  R  .  ■  1170, 

the  unstable  frequency,  at  l^/v  -  4.13  x  105(ra-1),  is  barely  outside  the 

unstable  region.  At  2vfe/U1  •»  0.20  and  R0  -  300,  the  unstable  frequency  is 

within  the  instability  region  for  the  mixing  layer  which  has  as  its  outer  bound 

2irf0/U1  ■>  0.27.  Comparison  of  the  unstable  frequencies  at  the  higher  Reynolds 

numbers  is  moot  since  they  are  regarded  as  behaving  nonllnearly.  However,  it  is 

evident  that  with  increasing  Reynolds  number  the  unstable  frequencies  expend 

further  into  the  damped  region  for  the  separated  layer.  As  for  the  mixing  layer 

comparison,  the  unstable  frequency  at  the  Reynolds  number  slightly  above 

critical  is  still  within  the  unstable  region,  but  at  l^/v  -  5.97  x  105(m-1)  and 
5  “1 

6.8l  x  10  (m  )  the  unstable  frequencies  are  outside  the  zone  of  instability. 

The  purpose  of  the  Toregoing  comparison  with  the  stability  calculations  is 
not  to  imply  any  direct  correspondence,  but  is  solely  for  the  purpose  of 
evaluating  whether  the  profiles  in  the  immediate  downstream  vicinity  of  the 
roughness  may  be  sufficiently  Inflectional  to  undergo  the  observed 
instabilities.  The  "strength”  of  the  inflection  and  lt3  position  relative  to 
the  surface  are  the  significant  criteria  in  assessing,  at  least  qualitatively, 
the  instability  of  an  inflectional  profile  in  a  wall-bounded  flow.  In  general, 
increased  proximity  of  the  inflection  to  the  surface  tends  to  have  a  stabilizing 
effect,  whereas,  an  increase  in  the  "strength"  of  the  inflection  i3 
destabilizing.  Inasmuch  as  it  has  been  previously  assumed  that  the  instability 
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la  essential  lnvlscld,  i.e.  that  the  Inflection  Is  at  the  critical  layer,  the 
Inflection  should  be  the  governing  criterion.  It  Is  seen  In  figures  43(a)  and 
(b)  that  the  measured  profiles  are  more  inflectional  than  either  the  Hartree  or 
mixing  layer  profiles,  and  It  can  be  concluded  that  the  measured  profiles  are 
sufficiently  inflectional  to  undergo  the  observed  instabilities.  In  this 
connection  it  would  be  helpful  to  have  a  stability  calculation,  albeit  two- 
dimensional,  of  a  profile  that  is  representative  of  the  measured  profile.  There 
is  also  an  indication  from  the  measured  profiles  in  figures  43(a)  and  (b)  that 
the  strength  of  the  inflection  may  increase  with  increasing  Reynolds  number. 

This  may  contribute  to  the  very  rapid  increase  in  amplification,  shown  in  figure 
42(a),  for  the  corresponding  range  of  Reynolds  numbers. 

The  foregoing  discussion  has  attempted  to  characterize  the  instability 
resulting  from  a  three-dimensional  roughness  element  as  an  inflectional 
Instability.  In  lieu  of  any  theoretical  calculations  for  evaluating  three- 
dimensional  instabilities  in  three-dimensional  flows,  the  rationalization  has 
taken  the  liberty  of  comparison  with  the  theoretical  guidelines  existing  for 
two-dimensional  instabilities  in  two-dimensional  flows,  and  a  reasonable  case 
for  inflectional  instability  has  been  made.  The  sensitivity  of  the  eddy 
shedding  frequencies  to  background  disturbances,  and  their  nondimenslonal 
scaling,  as  discussed  in  Section  3*2,  also  provide  Indirect  evidence  for  an 
inflectional  instability.  The  consideration  of  an  inflectional  instability 
provides  a  simple  model  for  the  behavior  of  a  three-dimensional  roughness 
element.  The  critical  roughness  Reynolds  number  reflects  the  separation  Into 
linear  and  nonlinear  instabilities.  In  the  linear  range,  at  subcritlcal 
Reynolds  numbers,  the  disturbance,  although  it  may  be  initially  amplified, 
cannot  sustain  Itself.  It  can  dissipate  laterally,  and  damp  as  it  travels 
downstream  into  more  stable  regions  where  the  mean  flow  profiles  become  less  and 
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leas  inflectional,  aa  evidenced  by  the  profiles  in  figure  41.  At  some  Reynolds 
number  the  instability  becomes  nonlinear.  In  this  range  the  disturbance  can  no 
longer  be  characterized  by  linear  theory.  A  frequency  which  may  be  stable 
according  to  linear  theory  may  well  be  unstable  when  it  is  nonlinear.  The 
amplification  in  the  nonlinear  range  of  Reynolds  numbers  Is  extremely  rapid,  and 
the  instability  results  in  a  rolled-up  vortex.  As  concluded  in  Section  3.3, 
this  vortex  is  intrinsic  to  the  transition  process.  The  final  onset  of 
turbulence  results  from  vortical  instability,  and  the  complex  vortex-vortex 
interaction  of  the  shedding  vortices  with  the  pre-existing  stationary  vortices. 

The  above  model,  however,  does  not  explicitly  account  for  the  observations 
of  Hall  (reference  17).  and  Morkovin  (reference  35)  referred  to  in  Section  3.3, 
i.e.  the  occurrence  of  "hairpin"  eddies  with  an  intensity  level  as  high  as  4 
percent  which  develop  and  decay  without  the  onset  of  turbulence.  This  behavior 
Implies  a  narrow  range  of  nonlinearity  with  Reynolds  number  for  which  the 
resulting  vortices,  and  vortex-vortex  interaction  are  not  unstable.  The 
intensity  level  of  4  percent  is  consistent  with  this  possibility,  and  the 
distributions  of  intensity  in  figure  42(a).  In  any  event,  this  aspect  merits 
further  investigation,  not  only  to  ascertain  its  role  In  the  context  of  the 
present  investigation,  but  in  view  of  its  importance  in  understanding  the  vortex 
dynamics  involved  in  the  onset  of  turbulence.  In  the  same  vein,  it  is  not 
evident  that  the  formation  of  the  "hairpin"  eddies  from  a  rolling  up  of  the 
spiral  vortex  filaments  can  be  reconciled  with  an  inflectional  instability.  If 
the  flow  visualization  is  to  be  taken  literally,  then  the  evolution  of  the 
"hairpin"  eddies  involves  a  vortical  instability.  On  the  other  hand,  their 
formation  from  an  inflectional  instability  would  Involve  nonlinear  development 
of  an  initial  Rayleigh  instability,  and  the  "hairpin"  eddies  may  well  be 
generated  by  a  time-dependent  concentration  of  spanwise  vorticity,  and  not  by  an 
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altered  state  of  the  spiral  vortex  filaments.  One  may  speculate  that  the 
vagaries  of  flow  visualization  do  not  provide  an  Insight  into  the  realities  of 
such  distinctions.  For  example,  is  the  flow  tracer  "captured"  by  one  set  of 
vortices  at  the  expense  of  the  other  thus  resulting  in  an  apparent  transition 
from  the  spiral  vortex  filaments  to  the  "hairpin"  eddies?  It  is  evident  that 
this  aspect  also  merits  further  investigation  in  order  to  arrive  at  a  proper 
understanding  of  the  induced  transition  to  turbulent  flow. 

The  critical  behavior  of  a  three-dimensional  roughness  element  in  inducing 
transition  to  turbulent  flow  is  clearly  demonstrated  by  the  distributions  of  U 
and  u'  at  the  various  Reynolds  numbers.  At  subcrltical  Heynolds  numbers  the 
distributions  return  to  that  which  would  exist  without  the  roughness  present. 

In  figure  44  the  mean  velocity  distributions  at  the  various  x-positions  ror  the 
subcrltical  Reynolds  numbers,  and  the  Reynolds  number  which  Is  slightly  above 
critical,  are  compared  nondimenslonally  with  the  mean  velocity  profile  measured 
without  the  roughness.  The  distance,  x,  required  to  return  to  the  undisturbed 
condition  is  dependent  on  the  Reynolds  number.  At  unit  Reynolds  numbers  of 
3.54  x  105,  and  4.13  x  105(m_1),  for  example,  the  recovery  occurs  in  x  -  12.7 
and  20.3  cm  respectively.  At  the  Reynolds  number  which  is  Just  above  critical 
the  mean  velocity  distributions  are  indicative  of  transitional  flow  consisting 
of  intermittent  turbulent  "spots".  The  changes  in  profile  shape  with  x  reflect 
an  increase  in  the  percent  of  time  that  the  flow  is  turbulent  with  increasing  x 
consistent  with  the  downstream  growth  of  turbulent  "spots".  However,  at  this 
Reynolds  number  the  flow  condition  is  extremely  sensitive  to  Reynolds  number. 
The  distributions  of  u*  at  this  Reynolds  number  reflect  this  sensitivity  in 
their  variation  with  x. 

As  discussed  in  Section  3.3  the  turbulent  "spots"  are  initiated  closer  to 
the  roughness  with  increasing  Reynolds  number  until  at  some  Reynolds  number  a 
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wedge  of  turbulence  may  be  regarded  as  being  attached  to  the  roughness.  Since  a 

finite  distance  is  required  for  the  instability,  the  wedge  of  turbulence,  albeit 

in  close  proximity  to  the  roughness,  Is  not  actually  attached.  At  a  unit 

Reynolds  number  of  6.81  x  1 0^(m  ^)t  which  is  the  highest  for  which  the 

distributions  were  measured,  the  corresponding  Rek  is  550,  and  the  wedge  of 

turbulence  has  its  origin  at  x  -  5.08  cm.  It  is  evident  from  the  progressive 

changes  in  the  distributions  of  U  and  u1  with  x,  at  the  two  highest  Reynolds 

numbers,  that  the  transition  to  a  turbulent  state  characteristic  of  a  "fully 

developed"  turbulent  boundary  layer  is  initiated  near  the  surface.  In  this 

connection,  it  should  be  noted  that  the  turbulence  at  the  origin  of  the  wedge, 

or  at  its  point  of  initiation  has  in  general  not  been  characterized.  In  the 

present  case,  turbulence  is  characterized  solely  by  the  u-f luctuation  having  a 

continuous  spectrum.  The  turbulence,  according  to  this  definition,  begins  near 

the  surface,  moves  out  in  y  with  increasing  x,  but  is  not  characteristic  of 

"fully  developed"  turbulence  until  x  -  30.5  cm.  This  is  seen  in  figure  45  in 

which  the  mean  velocity  distributions  at  various  x  for  unit  Reynolds  numbers  of 

5.97  x  10  and  6.81  x  10  (m  )  are  compared  with  a  turbulent  boundary' layer 

profile  measured  by  Purtell  et  al  (reference  42) .  The  position  at  which  the 

turbulence  in  the  wedge  becomes  characteristic  of  a  "fully  developed"  boundary 

layer  Is  examined  in  figure  46  in  which  the  distributions  of  u'  at  x  -  20.3  and 

5  *1 

30.5  on  at  a  unit  Reynolds  number  of  6.81  x  10  (m  )  are  compared  with  a 
distribution  from  reference  42  for  a  turbulent  boundary  layer  at  about  the  same 
Reynolds  number.  The  measurements  of  Purtell  et  al  were  made  with  a  linearized 
hot-wire  response  whereas  the  present  measurements  were  not.  However,  the 
magnitude  of  the  difference  between  the  distributions  is  such  that  attributing 
the  difference  to  the  method  of  measurement  is  questionable.  Either,  it  is  due 
to  other  experimental  uncertainties,  or  is  real.  The  latter  is  a  strong 
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possibility.  The  equivalence  of  Reynolds  numbers  based  In  one  ease  on  velocity, 
and  in  the  other  on  distance  has  not  been  thoroughly  evaluated  for  a  turbulent 
boundary.  Involved  herein  is  the  question  of  "new"  vs.  "aged"  turbulence.  The 
present  measurements  were  made  at  a  much  larger  U1  than  the  measurements  of 
Pur tell  et  al,  whereas  the  latter  were  made  at  what  may  be  considered  a  much 
larger  x.  In  view  of  these  uncertainties  it  is  concluded  that  the  turbulence  as 
well  as  the  mean  velocity,  are  characteristic  of  a  "fully  developed"  turbulent 
boundary  layer  at  x  -  30.5  cm.  The  values  of  x  are  not  scaled  with  k  since 
there  is  not  adequate  evidence  that  such  scaling  is  appropriate.  The  position 
in  the  wedge  at  which  the  turbulence  becomes  characteristic  of  a  "fully 
developed"  boundary  layer  is  observed  in  the  present  study  to  move  upstream  with 
increasing  Reynolds  number  but  whether  an  asymptotic  x/k  is  reached  has,  to 
date,  not  been  established  nor  has  the  effect  of  roughness  size,  shape,  and 
position  on  such  scaling  been  determined.  This  question  of  scaling  merits 
further  investigation,  and  it  is  apparent  that  considerable  care  should  be  given 
to  the  state  of  a  "tripped"  boundary  layer  In  order  to  ascertain  the  degree  to 
which  it  is  free  of  distortion. 

4.  Summary  and  Conclusions 

An  experimental  investigation  is  described  which  is  directed  toward 
extending  the  technical  database  and  Improving  our  understanding  of  the  manner 
in  which  single  three-dimensional  roughness  elements  Induce  earlier  transition 
from  laminar  to  turbulent  flow  in  boundary  layers.  The  Investigation  was 
carried  out  primarily  with  hemispherical  roughness  elements  in  a  well 
characterized  zero-pressure  gradient  laminar  boundary  layer  on  a  flat  plate. 

The  critical  roughness  Reynolds  number,  Re^,  at  which  turbulence  may  be  regarded 
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as  originating  from  the  roughness,  was  experimentally  determined  for  the  hemis¬ 
pherical  elements  herein  considered  to  be  Independent  of  roughness  size  and 
position.  The  value  of  Rek  based  on  the  roughness  height  and  the  velocity  in 
the  boundary  layer  at  the  height  of  the  roughness  was  325.  The  dependence  of 
Re^  on  roughness  shape  and  aspect  ratio,  as  well  as  its  imprecision  in  correlat¬ 
ing  the  behavior  of  single  roughness  elements,  is  demonstrated.  It  la  shown 
that  the  critical  behavior  of  three-dimensional  roughness  elements  and  the 
associated  Rek  are  a  consequence  of  the  critical  onset  of  an  eddy  generating 
process.  The  effect  of  a  steady  and  an  oscillatory  free-stream  velocity  on  the 
eddy  shedding  behavior  was  also  investigated.  The  results  demonstrate  that  the 
eddy  shedding  frequency  varies  linearly  with  the  square  of  the  Tree-stream 
velocity,  and  is  also  functionally  dependent  on  the  ratio  of  roughness  height  to 
position.  In  this  connection,  it  appears  that  the  boundary  layer  displacement 
thickness  at  the  roughness  position,  without  the  roughness  present,  is  a  more 
appropriate  length  scale  for  characterizing  the  Strouhal  behavior  than  the 
customarily  used  roughness  height.  It  Is  also  shown  that  roughness  shape  and 
aspect  ratio  influence  the  Strouhal  behavior.  It  is  concluded  that  the  eddies 
resulting  from  the  eddy  shedding  process  are  essential  and  directly  Involved  In 
the  transition  process.  They  may  be  regarded  as  being  elemental  "structures"  of 
the  turbulence  generated  by  the  complex  interaction  of  the  shedding  vortices 
with  the  pre-existing  stationary  vortices  that  lie  near  the  surface  and  are 
Inherent  to  the  flow  about  a  three-dimensional  obstacle  in  a  laminar  boundary 
layer.  The  mean  and  fluctuating  flow  fields  were  measured  on  the  centerline  at 
different  positions  downstream  of  a  hemispherical  roughness  element  which  had  a 
height  of  t.T  mm  and  was  positioned  at  91.4  cm  from  the  leading  edge  of  the  flat 
plate.  The  measurements  were  made  at  Reynolds  numbers  which  were  below  and 
above  critical.  Such  measurements  not  only  extend  the  technical  database  for 
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characterizing  roughness  element  behavior,  but  also  illustrate  that  considerable 
care  should  be  given  to  the  state  of  the  "tripped”  boundary  layer  in  order  to 
ascertain  the  degree  to  which  it  is  free  of  "distortion".  These  measurements 
also  permitted  within  the  limitations  Imposed  by  present  stability  theory,  an 
evaluation  of  whether  an  inflectional  instability  may  be  the  fundamental 
mechanism  underlying  the  behavior  of  three-dimensional  roughness  elements.  It 
is  concluded  that  the  behavior  of  three-dimensional  roughness  elements  is 
consistent  with  such  an  instability.  However,  further  investigation  is  required 
in  order  to  reconcile  this  view  with  the  observations  from  the  flow 
visualization  studies. 
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c.  Primary  and  secondary  duets 
Figure  1.  Schematic  drawing  of  Unsteady  Flow  Facility. 
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Figure  Photographs  of  the  overall  experimental  arrangement. 
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a.  Traversing  mechanism 


b.  Roughness  element,  and  hot-wire  probe 
Figure  5.  Close-up  view  of  experimental  arrangement. 
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Figure  8.  Intensity  of  u-fluctuation  without  roughness. 
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Figure  10.  Spectrum  of  u-fluctuation  in  the  freestream 
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Figure  11.  Spectrum  of  the  output  of  an  accelerometer  attached  to  the  fan 
housing. 


AEDC  TR-87-7 


Figure  12.  Variation  of  Reynolds  number  of  transition  with  roughness 
Reynolds  number  for  hemispherical  roughness  elements. 
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Figure  14.  Oscillograms  of  u-fluctuation  downstream  of  a  hemispherical 
roughness  element  illustrating  variation  in  frequency  wj.th 
freestream  velocity;  k  =  1.7  mm,  x|<  =  91.4  cm,  x  =  2.54  cm,  z  = 
0.0  cm. 
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U,  =  7.0  m/s 
y,  =  1  96  mm 
y2  =  3.86  mm 


11,=  8.2  m/s 
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U  ,=  10. 1  m/s 
y,  =  2.79  mm 
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y,  =  2.41  mm 
y2—  4  52  mm 


U  ,=  8.3  m/s 
y,  =  2.92  mm 
y2=  5.03  mm 


^  rf'rrrrrr^rrrf'\ 


U  ,=  10  2  m/s 
y,  =  3.76  mm 
y2  =  5.87  mm 


U  ,=  1  1.7  m/s 
y,  =  3. 18  mm 
y2  =  5.28  mm 


Figure  15.  Oscillograms  of  simultaneous  u-fluctuations  downstream  of  a 
hemispherical  roughness  element  illustrating  phase  in  y- 
direction;  k  =  1.7  mm,  x^  =  91. ^  cm,  z  =  0.0  cm. 
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Figure  16.  Oscillograms  of  simultaneous  u-fluctuation  downstream  of  a 
hemispherical  roughness  element  illustrating  phase  in  z- 
direction;  k  =  1.7  mm,  x^  =  1.27  cm,  x  =  1.27  cm,  y  =  3.18  ran. 
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Figure  17-  Variation  of  eddy  shedding  frequency  with  freestrean  velocity  for 
hemispherical  roughness  elements  of  different  size  and  position. 
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Figure  19-  Strouhal  behavior  for  hemispherical  roughness  elements  with 
boundary  layer  scaling. 
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Figure  20.  Variation  of  k/8*k  with  roughness  Reynolds  number  for 
hemispherical  roughness  elements  of  different  size  and  position. 


Figure  22.  Comparison  of  variation  of  eddy  shedding  frequency  with 
freestream  velocity  for  a  hemispherical  and  a  cylindrical 
roughness  element. 
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Figure  25.  Comparison  of  Strouhal  behvaior  for  a  cylindrical  roughness 
element  observed  in  the  present  study  with  the  observations  of 
Strouhal  behavior  in  reference  18. 
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Figure  26.  Phase  and  amplitude  response  of  the  laminar^  boundary  layer 
without  roughness  to  a  frees tr earn  oscillation;  Ui  =  6.8  m/s,  x.  = 
9^  cm,  n  =  1.0  Hz. 
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U,  =  7.6  m/s  AU,  =  0.98  m/s  n  =  1.0  Hi  1 

Figure  28.  Strip-chart  recordings  of  simultaneous  u-fluctuations  in  an 
oscillatory  freestream,  and  in  the  boundary  layer  downstream  of  a 
hemispherical  roughness  element;  k  =  1.7  nit,  xk  =  91.4  cm,  x  = 
2.54  cm,  y  =  3.18  mm,  z  =  0.0  cm. 
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0,  =  5.8  m/s  AU,  =  0.42  m/s  n  =  2.0  Hz 


U,  =  7.6  m/s  AU,  =  0.64  m/s  n  =  2.0  Hz 

Figure  29.  Strip-chart  recordings  of  simultaneous  u-fluctuations  in  an 
oscillatory  freestream,  and  in  the  boundary  layer  downstream  of  a 
hemispherical  roughness  element;  k  =  1.7  mm,  xfc  =  91.4  cm,  x  = 
2.54  cm,  y  =  3.18  mn,  z  =  0.0  cm. 
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Figure  29 
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0.89  m/s  n  =  2.0  Hz 
Concluded. 
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U,  =  8.3  m/s  AU,  =  0.00  m/s  n  =  0.0  Hz 


0,  =  8.1  m/s  AU,  =  1.08  m/s  n  =  1.0  Hz 


0,  =  9.4  m/s  AU,  =  1.34  m/s  n  =  1.0  Hz 

Figure  32.  Strip-chart  recordings  of  simultaneous  u-fluctuations  in  an 
oscillatory  freestream,  and  in  the  boundary  layer  downstream  of  a 
hemispherical  roughness  element;  k  =  1.7  mm,  x^  =  91.4  cm,  x  = 
30.5  cm,  y  =  6.35  mm,  z  =  0.0  cm. 
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0,  =  7.5  m/s  All,  =  0.63  m/s  n  =  2.0  Hz 
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U,  =  8.2  m/s  AU,  =  0.00  m/s  n  =  0.0  Hz 


0,  =  7.2  m/s  AU,  =  0.90  m/s  n  =  1.0  Hz 


0,  =  8.2  m/s  AU,  =  1.10  m/s  n  =  1.0  Hz 


D,  =  9.6  m/s  AU,  =  1.36  m/s  n  =  1.0  Hz 

Figure  34.  Strip-chart  recordings  of  simultaneous  u-fluctuations  in  an 
oscillatory  freestream,  and  in  the  boundary  layer  downstream  of  a 
hemispherical  roughness  element;  k  =  1.7  nm,  xi<  =  91.^  cm,  x  = 
61.0  cm,  y  =  0.61  mm,  z  =  0.0  cm. 
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U,  =  9.4  m/s  AU,  =  1.32  m/s  n  =  1.0  Hz 

Figure  36.  Strip-chart  recordings  of  simultaneous  u-fluctuations  in  an 
oscillatory  freestream,  and  in  the  boundary  layer  downstream  of  a 
hemispherical  roughness  element;  k  =  1.7  mn,  *k  =  91.4  cm,  x  = 
61.0  cm,  y  =  6.35  mm,  z  =  0.0  cm. 
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0,  =  9.4  m/s  AU,  =  0.88  m/s  n  =  2.0  Hz 

Figure  37.  Strip-chart  recordings  of  simultaneous  u-fluctuations  in  an 
oscillatory  freestream,  and  in  the  boundary  layer  downstream  of  a 
hemispherical  roughness  element;  k  =  1.7  mm,  =  91.4  cm,  x  = 
61.0  cm,  y  =  6.35  mm,  z  =  0.0  cm. 
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Figure  38.  Comparison  of  the  duration  of  eddy  shedding  with  the  duration  of 
resulting  turbulent  "spot." 
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Figure  40.  Spanwise  distributions  of  the  intensity  of  u-fluctuation  at 
various  positions  f ran  the  surface  for  varying  unit  Reynolds 
number;  k  =  1.7  nm,  x|i  =  91.4  cm,  x  =  2.54  cm. 
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y  (mm) 

b.  x  =  2.54  cm 
Figure  41.  Continued. 
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g.  x  =  30.5  cm 
Figure  41.  Concluded. 
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a.  7  =  1.27  cm 

Figure  M2.  Intensity  of  u-fluctuation  for  varying  unit  Reynolds  number;  k  = 
1.7  nm,  xk  =  91-4  cm,  z  =  0.0  cm. 
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Figure  42(b) 


y  (am) 

b.  x  =  2.54  cm 
Figure  42.  Continued 
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d.  x  =  7.62  can 
Figure  42.  Continued. 
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b.  Separated  boundary  layer 

Figure  43.  Comparison  of  mean- velocity  profiles  at  1,27  cm  downstream  of  a 
hemispherical  roughness  element  with  the  mean-velocity  profiles; 
k  =  1.7  mm,  *k  =  91.4  cm,  z  =  0.0  cm. 
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Figure  46.  Comparison  of  intensity  of  u-fluetuation  downstream  of  a 
hemispherical  roughness  element  with  the  intensity  of  u- 
fluctuation  for  a  "fully-developed"  turbulent  boundary  layer;  k  = 
1.7  nm,  xfc  =  91.4  cm,  z  =  0.0  cm,  Ui/v  =  6.81  x  105  (m  ). 
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